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INTRODUCTION 

Utilization  of  electricity  on  the  farm  and  in  the  farm  home  is  increas- 
ing along  with  the  development  of  new  equipment  as  a  result  of  further 
application  of  basic  principles  which  were  discovered  50  to  100  or  more 
years  ago.   One  such  electro-physical  principle,  known  as  thermoelectricity, 
has  created  interest  in  recent  years  with  regard  to  the  use  of  the  Peltier 
effect  for  cooling  and/or  heating  purposes.   This  is  the  reverse  process  to 
thermoelectric  power  generation  (Seebeck  potential)  of  direct  current  by 
heating  alternate  junctions  of  pairs  of  dissimilar  metals  or  semi-conductors, 
the  latter  being  required  for  high  efficiency.   Utilization  of  the  refriger- 
ating affect  achieved  in  this  manner  has  been  limited  to  laboratory  experi- 
ments due  to  the  minuscule  nature  of  the  Peltier  effect. 

A  French  physicist,  Jean  C.  A.  Peltier  in  1834  first  observed  the 
reverse  of  the  thermoelectric  phenomenon  discovered  in  1832  by  Thomas  J. 
Seebeck,  a  German  physicist.   At  that  time  it  would  have  been  theoretically 
possible  to  generate  electricity  with  thermocouples  at  an  efficiency  of  about 
3  percent,  which  would  have  been  as  good  as  the  best  steam  engines  available. 

During  the  past  decade  or  so  progress  in  the  development  of  thermoelec- 
tric materials  have  produced  more  efficient  semi-conductor  compounds.   Lead 
Telluride  is  one  such  compound  that  is  now  being  used  in  such  applications 
as  water  coolers,  ice  cube  makers,  air  conditioners,  and  numerous  instrumen- 
tation devices.   The  material  development  has  reached  a  plateau  during  the 
past  several  years.   However  the  method  of  assembling  these  materials  into 
thermoelectric  junctions  has  continued  to  advance  and  thus  further  improve 
the  efficiency  of  the  devices. 

There  is  also  continued  advancement  in  mass  production  techniques  for 


assembling  various  numbers  of  thermoelectric  junctions  into  a  module.   Such 
techniques  will  further  reduce  the  cost  per  Btu  of  cooling  capacity  which 
will  make  larger  thermoelectric  devices  more  competitive  with  conventional 
refrigeration  systems. 

The  reverse  cycle  refrigeration  system  or  "heat  pump"  as  it  is  commonly 
called,  is  capable  of  transferring  heat  from  a  low  energy  level  to  a  high 
energy  level.   This  device  can  transfer  more  heat  (Btu/Hr)  than  the  heat 
equivalent  (Btu/hr)  of  the  system  power  input.   The  ratio  of  the  heat  trans- 
ferred to  the  heat  equivalent  of  the  power  input  is  called  the  coefficient 
of  performance  or  C.O.P.  of  the  device.   A  thermoelectric  "heat  pump"  invar- 
iably has  a  heating  coefficient  of  performance  greater  than  unity.  At  the 
present  state  of  thermoelectric  material  development,  the  cooling  coefficient 
of  performance  is  usually  less  than  unity.   However  it  is  predicted  that  as 
more  efficient  thermoelectric  materials  are  developed  and  better  design 
methods  are  found  for  applying  them  in  heat  transfer  devices,  a  cooling 
coefficient  of  performance  greater  than  unity  can  be  accomplished. 

Also  this  transfer  of  heat  can  be  accomplished  by  thermoelectric 
devices  which  have  no  moving  parts.   The  refrigeration  effect  and  heat 
transfer  is  accomplished  solely  by  the  flow  of  electrons  through  the  thermo- 
electric materials.   Today's  agricultural  scientists  apply  new  principles 
and  methods  of  heat  transfer  to  their  design  needs.   The  improvement  of 
conventional  systems  by  eliminating  moving  parts  must  be  considered  to  be 
of  major  importance.   Several  theoretically  possible  thermoelectric  appli- 
cations in  agriculture  are  cooling  milk  and  heating  water  (or  air);  milk 
cooling  and  pasteurization;  humidity  and  dew  point  controls;  equipment 
operator  air  conditioning;  animal  body  cooling  or  heating;  dehumidif iers; 


electronic  control  element  cooling,  water  cooling;  refrigeration  of  semen, 
antibiotics  and  vaccines  and  small  isolated  power  supplies. 

THE  PROBLEM 

The  design  of  thermoelectric  heating/cooling  systems  has  heretofore 
been  investigated  for  systems  that  generally  have  somewhat  fixed  conditions 
on  both  sides  of  the  heat  transfer  device.   Such  systems  are  represented  by 
water  coolers,  refrigerators,  air  conditioners,  ice  cube  machines  and 
freezers. 

Methods  of  estimating  the  cooling  capacity,  power  requirements  and 
other  factors  have  utilized  complicated  charts  and  computer  programs.   There 
appears  to  be  very  little  information  on  engineering  design  equations  for 
determining  the  relationship  of  the  important  variables  in  a  transient 
system.   Such  a  device  can  be  represented  by  a  milk  cooler/water  heater  or 
a  milk  pasteurizer/water  cooler,  milk  cooler/water  heater  combination. 

It  would  appear  that  a  general  engineering  design  equation  which 
relates  the  important  variables  would  be  a  useful  tool  for  the  engineer. 

OBJECTIVES 

The  objectives  of  the  study  were: 

1.  To  investigate  the  development  of  a  general  prediction  equa- 
tion that  describes  the  relationship  of  the  important  vari- 
ables in  a  thermoelectric  heating/cooling  system. 

2.  To  verify  the  general  prediction  equation  by  laboratory 
experiments. 

The  experimental  research  was  planned  to  incorporate  two  flat-plate 


liquid  heat  exchangers  with  a  12-couple  thermoelectric  module  mounted 
between  these  heat  exchangers.   The  thermoelectric  module  is  of  a  type 
that  is  commercially  available  at  this  time. 

REVIEW  OF  LITERATURE 

Only  one  reference  has  been  found  regarding  the  use  of  thermoelectric- 
ity in  agriculture.   Golubyatnikov  (10)  reports  on  a  thermoelectric  milk 
cooler/water  heater  unit  being  tested  in  Russia.   He  states, 

Its  relevant  data  are:   throughput  100-150  lit./hr.  (26-40 
gal./hr.)  (or  approximately  220-340  lbs./hr.)  with  milk  storage 
capacity  up  to  2000  Kg  (510  gal.);  refrigerating  capacity,  2700- 
3000  Kcal./hr.  (10,700-11,900  Btu./hr.);  cooling  coefficient  E 
(CO. P.)  =  1.36  -  1.92;  power  requirement  =  0.9-2.7  KW,  d.c. 
current  I  =  140  -  150  A  (amperes);  milk  can  be  cooled  to  a  tempera- 
ture of  8-5°C  (46-4l°F) ;  the  energetic  characteristic  of  the  semi- 
conductor material  Z  =  1.70  x  10""-VoC;  the  weight  of  semiconductor 
material  used  is  6.0-7.5  Kg.  (13-16.5  lbs.);  the  number  of  thermo- 
couples:  up  to  1000. 

The  thermoelectric  refrigerator  constitutes  a  thermoelectric 
battery  of  circular  high-current  thermocouples  consisting  of  a 
ternary  alloy  Te  Bi  Se  negative  branch  and  a  Te  Bi  Sb  positive 
branch,  connected  in  series.   It  was  found  that  the  inertia  of 
circular  thermocouples  is  small,  and  that  they  can  provide  a  high 
power  output  with  stable  characteristic. 

By  this  circular  thermocouple  design,  resulting  in  much 
simplified  manufacturing  technology,  semiconductor  power  systems 
can  be  realized  .  .  . 

Golubyatnikov  further  describes  the  thermoelectric  rectifier  power  supply 

water  cooler.   Thermostats  control  the  milk  input  and  output  for  the  device. 

The  milk  is  first  pasteurized  and  then  cooled.   The  heat  transfer  mechanism 

used  provides  a  solution  of  the  problem  of  the  cooling  of  milk  with  material 

having  a  small  figure  of  merit. 

He  states  further, 

The  counter-convective  heat  exchange  is  effected  (provided) 
by  the  appropriate  direction  of  the  milk  and  cooling  water  flows, 


and  by  passing  the  liquid  to  be  cooled  over  the  radiator  surface 
of  the  generator  (thermoelectric  cooler).   This  assures  maximum 
heat  loss  and  a  maximum  heat  transfer  coefficient  and  consequently 
an  increase  in  generator  efficiency.   Thus  an  energetic  characteris- 
tic Z  =  1.70  x  10"3/°C  obtained  by  the  USSR  Semiconductor  Institute, 
is  perfectly  adequate  for  solving  the  thermal-engineering  and 
hydromechanical  problems  connected  with  the  design  of  an  economi- 
cally operating  thermoelectric  refrigerator-heater  for  dairy  farms. 

A  paper  prepared  by  Charles  L.  Feldman  et  al.  (7)  explains  a  method  by 

which  cascaded  peltier  coolers  can  be  optimized  under  certain  specified 

conditions  such  as  hot  junction  temperature.   The  figure  of  merit  was  assumed 
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to  be  2.8  x  10   K  '  and  T,  =  300K.   A  low  temperature  freezer  was  designed 

from  theoretical  calculations  using  bismuth  telluride  thermoelectric  modules. 

A  report  (1)  by  Dr.  E.  F.  Cox  of  the  Whirlpool  Corporation  states, 
"Except  for  economics,  thermoelectric  heat  pumping  should  be  excellently 
suited  to  distillation  processes."   He  presents  an  abbreviated  and  simplified 
theoretical  analysis  to  determine  the  optimum  thermoelectric  material  pro- 
perties to  maximize  efficiency  of  salt  water  still  operation.   Several  ideas 
are  presented  for  making  short  thermoelements  and  junctioning  them  to 
bridges.   The  possibility  of  thin  films  of  metals  or  semiconductors  are 
discussed  and  studies  along  these  lines  are  recommended  as  a  wide-open  area 
for  experimental  research. 

A  series  of  reports  by  the  United  States  Naval  Research  Laboratory 
Washington,  D.C.  (2)  review  the  status  of  thermoelectric  materials  and 
device  developments.   The  final  report  (NRL  Memorandum  Report  1404)  in 
this  series  discusses  briefly  Peltier  heat  pumps,  parameters  and  device 
possibilities.   The  report  states,  "Materials  developments  now  make  pos- 
sible the  construction  of  practical  reversible  Peltier  heat  pumps  and  simple 
refrigerators  for  a  variety  of  military  applications." 

D.  C.  Siegla  et  al.  (19)  have  analyzed  the  heat  transfer  characteristics 


of  a  thermoelectric  module  for  refrigeration.   A  one  dimensional  and  a  three 
dimensional  analysis  are  described.   In  both  analyses  the  relationship  of 
junction  temperature  and  dimensions  of  the  thermoelectric  element  are 
evaluated  with  respect  to  the  heat  transfer  capability  of  the  module. 

In  several  progress  reports  by  the  Carrier  Corporation  to  the  Depart- 
ment of  the  Navy  (22),  research  is  described  relating  to  design  of  submarine 
air  conditioning  units.   An  analog  computer  study  is  described  which  program 
produces  operation  data  based  on  selected  values  of  junction  temperature 
differences.   Thermoelectric  submarine  air  cooling  systems  utilizing  sea 
water  exchangers  are  estimated  to  be  71  percent  as  heavy  as  a  conventional 
reciprocating  refrigeration  system. 

A  NRL  status  report  by  Dr.  J.  W.  Davisson  and  Joseph  Pasternak  (3) 
describes  theoretical  analysis  by  analog  computer  programs  of  thermoelectric 
devices  for  use  in  submarines.   In  each  design  certain  factors  were  held 
constant  such  as  sink  water  temperature,  air  flow  through  cooling  coil  or 
cooling  capacity.   Therefore  all  the  important  variables  are  not  combined 
into  one  solution. 

Another  NRL  status  report  by  Dr.  J.  W.  Davisson  and  Joseph  Pasternak 
(4)  summarizes  thermoelectric  materials  characteristics  and  includes  such 
parameters  as  figure  of  merit,  Seebeck  coefficient,  electrical  resistivity, 
thermal  conductivity  and  others. 

A  report  by  E.  W.  Frantti,  Westinghouse  Electric  Corporation  (8) 
describes  the  analysis  of  a  thermoelectric  air  conditioning  system  for 
submarines.   Test  loops  were  set  up  to  measure  the  operating  characteris- 
tics of  a  thermoelectric  module  in  the  laboratory.   Various  relationships 
were  derived  from  laboratory  measurements.   However  the  relationships  were 


not  of  a  mathematical  form  such  that  the  variables  could  be  combined.   For 
certain  fixed  conditions,  curves  are  presented  which,  for  example,  show  the 
relationship  of  heat  pumping  rate  and  CO. P.,  to  current,  chill  water  flow 
or  heat  sink  water  flow. 

The  transient  performance  of  a  thermoelectric  couple  under  step-current 
control  is  described  in  a  paper  by  W.  F.  Stoecker  and  J.  B.  Chaddock  (20). 
The  temperature  of  the  cold  junction  with  respect  to  time  is  shown  following 
a  step-change  in  current.   For  example  this  characteristic  is  shown  for 
step-change  from  0  to  10  amperes  and  from  10  to  20  amperes  and  vice  versa. 
In  all  cases  the  hot  junction  temperature  was  held  constant  at  310  K. 

Magazine  articles  by  Alwin  B.  Newton  (15)  discuss  the  selection  of 
materials  and  prediction  of  performance  for  thermoelectric  systems.   Conven- 
tional charts  are  presented  which  show  performance  for  a  fixed  sink  tempera- 
ture.  Analysis  of  other  heat  transfer  variables  is  not  included.   The 
author  points  out  a  recent  development  of  ceramic  interfaces  between  the 
module  junctions  and  the  surface  being  cooled  or  heated.   This  development 
has  reduced  the  temperature  drop  across  such  interfaces  from  about  20  to 
25  degrees  F.  to  about  an  average  of  2  F.   Power  and  control  systems  are 
also  discussed. 

Rajic  presents  in  his  article  Refrigeration  Techniques  and  Thermo- 
electric Phenomena  (18),  the  theoretical  bases  of  analyzing  thermo-electric 
phenomena,  and  of  the  technical  calculations  of  thermo-electric  refrigera- 
tion devices.   To  determine  characteristics  of  thermo-electric  elements, 
scientists  have  suggested  various  coefficients,  such  as:   "effective  thermal 
force"  (E.  Altenkirch),  dimensionless  coefficient  0  (H.  J.  Goldsmid,  R.  W. 
Douglas)  and  "effectiveness  of  thermoelement"  Z  (Ioffe).   Z=  <*   ©  /K 


where  <?C  is  the  coefficient  of  proportionality  of  the  thermo-electric  excit- 
ing force  according  to  Zebeke;  E  =  OC  (T  -  T  )  where  E  is  the  voltage,  T  and 
T  -  the  temperature  of  the  hot  and  cold  junction  respectively;  £    -  the 
specific  electric  conductivity  of  the  thermo-electric  element,  and  }\    is  the 
coefficient  of  its  heat  conductivity.   If  these  characteristics  are  known, 
it  is  easy  to  determine  the  maximum  value  of  the  hot  and  cold  junction 
temperature  (T/T  )    for  obtaining  the  maximum  value  of  operation  current 

in  the  circuit,  I      .   The  efficiency  of  the  refrigerating  unit  is 
optim.  J  o      & 

determined  by  the  refrigeration  coefficient  £  found  from  equation 

w  max 

£max  =  T/(T-To)  =  E  VI  +  Z(T+TQ)/2  -  T/T  ]/|yi  +  Z(T+TQ)/2  +  l].   Data  are 
given  on  the  selection  of  materials  for  thermo-electric  elements,  in  par- 
ticular, some  characteristics  of  semiconductor  alloys  Bi  Te   (52.8%  Bi  + 
47.8%  Te)  and  Sb  Te,  (38.9%  Sb  +  61.1%  Te)  which  were  studied  at  the 
Institut  Poluprovodnikov  AN  USSR  (Institute  of  Semiconductors,  AS  USSR). 
The  conclusions  indicate  advantages  of  thermo-electric  refrigeration  over 
the  classical  liquid  one.   For  further  successful  utilization  of  thermo- 
electric phenomena  in  refrigeration  engineering,  joint  efforts  of  physicists, 
chemists  and  metallurgists  are  imperative. 

A  paper  by  R.  C.  Miller  et  al.  (14)  theoretically  analyzes  the  depend- 
ence of  the  figure  of  merit  (Z)  of  thermoelectric  material  on  the  energy 
bandwidth.   It  was  found  that  Z  approaches  zero  as  the  bandwidth  goes  to 
zero.   The  existence  of  an  optimin  bandwidth  is  established.   However  the 
magnitude  could  not  be  ascertained. 

A  paper  by  Roland  W.  Ure,  Jr.  (21)  presents  comments  on  thermal  con- 
ductivity in  two-phase  alloys.   Quantitative  agreement  between  the  predic- 
tions of  the  model  used  and  experimental  data  were  obtained. 


A  paper  by  I.  N.  Pomazanov  and  P.  L.  Tikhomikov  (Russia)  (16)  discusses 
the  operation  of  a  thermoelectric/generator/ refrigerator  operating  from  a 
gas  burner.   With  a  power  output  of  from  2  to  3KW,  sufficient  refrigeration 
capacity  was  obtained  to  freeze  about  0.7  kg  (1.54  lbs.)  of  water  per  hour. 
One  operational  disadvantage  of  the  system  described  was  the  high  rate  of 
flow  of  cooling  water  (100-120  liters/hr.)  (26.4-31.7  gal./hr.). 

In  a  report  to  the  Department  of  the  Navy  (23)  Carrier  Corporation 
research  work  on  thermoelectric  air  coolers  and  water  chillers  are  described, 
Exact  differential  equations  were  written  for  heat  and  mass  transfer  rela- 
tionships and  these  equations  along  with  basic  thermoelectric  equations  were 
solved  by  an  analog  computer. 

In  these  analyses  fixed  cooling  compartment  temperatures  were  assumed. 
Figure  of  merit  of  thermoelectric  materials  ranged  from  2.75  x  10  /  C  to 
3.04  x  10  /  C.  An  air  cooler  of  about  12,000  Btu/hr.  capacity  was 
designed  and  tested. 

A  report  by  the  Whirlpool  Corporation  (24)  to  the  Department  of  the 
Navy  describes  the  design  of  a  specific  prototype,  including  controls,  of 
a  thermoelectric  refrigerating  system  for  use  aboard  submarines.   One  room 
is  held  at  0  F  and  another  can  be  held  at  33  F  or  0  F  as  may  be  required. 
It  was  assumed  that  thermoelectric  materials  would  have  a  figure  of  merit 
of  at  least  2.5  x  10  /  C.   It  was  recommended  from  this  study  that  a  proto- 
type refrigeration  system  be  assembled  and  tested. 

A  report  by  the  York  Corporation  to  the  Department  of  the  Navy  (25) 
includes  a  thorough  study  of  present  food  preservation  and  storage  practices. 
New  food  preservation  methods  such  as  freeze  dehydration  are  discussed  and 
revised  food  storage  tables  were  developed.   Another  part  of  this  report 
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presents  a  preliminary  study  of  the  application  of  thermoelectricity  to 
submarine  food  storage  facilities. 

Comments  by  W.  Maurice  Pritchard  in  the  Proceedings  of  the  IEEE  (17), 
refer  to  previous  work  on  calculating  and  optimizing  the  coefficient  of 
performance  of  thermoelectric  cooling  devices.   Previous  studies  have  been 
based  on  the  average  thermoelectric  properties  and  have  neglected  the 
Thomson  effect.   Also  the  effects  of  junction  resistance  and  a.c.  ripple 
in  the  power  supply  are  important  factors  to  be  considered.   The  writer 
presents  a  mathematical  approach  to  the  optimization  problem. 

A  status  report  by  Dr.  J.  W.  Davisson  and  Joseph  Pasternak  (5)  presents 
thermoelectric  materials  data.  The  authors  point  out  that  no  dramatic  gains 
in  the  figure  of  merit  have  been  accomplished  at  the  time  of  their  report. 

However  they  state  that  Bell  Laboratories  have  developed  a  bismuth  antimony 

-3 
alloy  with  a  Z  of  better  than  5  x  10   at  80  K.   This  is  a  very  high  Z  and 

is  a  usable  performance  at  cryogenic  temperatures. 

A  paper  by  Marvin  A.  Fuller  (9)  discusses  the  mass  production  tech- 
niques of  thermoelectric  modules.   He  points  out  that  a  simple  flat  con- 
figuration which  can  be  sandwiched  between  two  heat  sinks  is  desirable. 
Junction  resistance  problems  are  discussed  along  with  the  process  apability 
level. 

A  report  by  Whirlpool  Corporation  to  the  Department  of  the  Navy  (26) 
describes  a  design  study  to  optimize  the  application  of  a  thermoelectric 
device  for  refrigeration  aboard  submarines.   The  study  clearly  showed  that 
it  is  possible  and  practical  to  thermoelectrically  provide  such  refrigeration 
in  food  storage  rooms. 
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THEORETICAL  ANALYSIS 
Dimensional  Analysis 

The  variables  which  influence  a  physical  system  can  often  be  correlated 
by  methods  of  dimensional  analysis.   Such  an  analysis  offers  a  means  of 
simplification  by  combining  these  factors  into  dimensionless  groups  which 
can  facilitate  experimental  research. 

Dimensional  analysis  establishes  the  relationship  of  these  variables  in 
a  model  which  can  then  be  compared  with  a  prototype  device. 

To  utilize  dimensional  analysis  the  variables  considered  to  be  impor- 
tant to  the  physical  system  are  first  selected.   Then  these  variables  are 
grouped  into  dimensionless  "77""  terms. 

The  "77""  terms  can  then  be  related  by  an  equation  of  the  form: 


irl  =  f(72,  wr...  ^n) 


In  order  to  determine  the  general  equation  the  nature  of  the  function 
"f"  must  be  established.   Satisfactory  predictions  can  be  made  if  the 
theoretical  model  and  prototype  are  operated  so  that  the  values  of  the  "77"" 
terms  in  each  case  are  equal. 

The  method  of  evaluating  the  function  "f"  is  to  hold  one  of  the  inde- 
pendent "77""  terras  constant  and  vary  the  other.   This  method  is  then  applied 
to  each  of  the  "77""  terms  and  the  resulting  component  equations  can  be 
combined  to  give  the  general  relationship.   The  final  equation  may  be  a 
polynomial  type  or  exponential  type  depending  on  the  relationship  of  the 
"77""  terms. 
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Dimensionless  Groups  and  "Jfu   Terms 

The  design  of  this  investigation  was  developed  around  the  principles  of 
dimensional  analysis.   The  formation  of  dimensionless  groups  was  investi- 
gated with  an  IBM  1620  computer  to  provide  theoretical  data.   This  data  was 
supplied  to  an  IBM  1410  computer  which  was  programmed  to  statistically 
analyze  the  data.   The  output  from  the  1410  computer  yielded  a  general 
prediction  equation  which  could  be  experimentally  investigated. 

In  analyzing  the  dimensionless  "77""  groups  the  following  variables  were 
considered  to  be  pertinent: 

Factor        Description  Dimension 

Mass-Length-Time 
Temperature  Units 


wb 
Ati 

At 

c 

W 
c 

W, 


Cooling  capacity 

Power  required 

Difference  between  inlet  water 
temperatures 

Drop  in  water  temperature  due 
to  cooling 

Water  flow  rate  through  cooling  sink 

Water  flow  rate  through  heating  sink 


2  -3 
ML  T  ° 

2  -3 

ML  T 


9 


0 


MT 


-1 


MT 


-1 


The  required  number  of  U7T"   terms  can  be  determined  as  follows: 


P   =  K(Q  )a,  <At.)\  (At  )C,  (W  )d,  (W,  )e 
wb      x:       l       c      c      h 


(1) 


Substituting  the  dimensions  for  each  term: 


ML 


ML 


2\a 


• 6".  •'•  $•  (ff 


(2) 
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Equations  can  then  be  written  by  summarizing  the  exponents  for  each 
dimension  as  follows: 

£M:  1  =  a  +  d  +  e  (3) 

2L:  2  =*  2a  (4) 

£T:  -3  =  -3a  -  d  -  e  (5) 

£3:  0  =  b  +  c  (6) 

From  these  equations  the  following  is  obtained: 

a  =  1 
d  +  e  =  0  or  d=-e 

c  =  -b  or  b  =  -c 
Substitution  of  the  above  values  into  equation  (2)  yields: 

From  equation  (7)  by  rearranging  terms  the  following  results: 

h  -  #  fey 

The  relationship  between  the  three  "77""  terms  thus  formed  can  be 
analyzed  by  holding  one  constant  while  the  other  two  are  varied.  The 
analysis  was  developed  as  follows: 

First  set  of  data  from  computer: 


"lm   f(^"2*   3*   wnere;   (77"  constant 

Q         At.        W 

77-  =-c-  7T  =  — i   7T=-£- 
'  1   P   '  "  2   At  '  "3   w\ 
wb         c         h 
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Second  set  of  data  from  computer: 

TTl   =  f (  7T   ,  T2)  (  7T2   constant) 

In  the  analysis  another  "jf"   term  was  used  as  a  level  indicator  to  cor- 
relate differences  in  cooling  capacity  and  water  inlet  temperatures  on  the 
cooling  side.   The  fourth  "77""  terra  was  as  follows: 

T 
77"  -  wic  _  inlet  water  temperature  on  cooling  side 

4    t     drop  in  water  temperature  due  to  cooling 

This  "77""  term  was  held  constant  at  four  different  values  of  inlet 
water  temperature  while  At  was  held  at  a  constant  value  of  2  F.   The 
relationship  between  the  other  "7f"   terms  was  then  determined  at  these 
levels.   The  data  were  then  combined  to  yield  a  general  prediction  equation. 
A  family  of  equations  of  surfaces  can  be  obtained  for  various  values  of  At  . 
Engineering  design  can  begin  by  selecting  a  value  of  At  and  then  using  the 
general  prediction  equation  the  other  pertinent  variables  can  be  related. 


Calculations  to  Determine  Theoretical  Equation  for 
Heat  Transfer  Across  the  Heat  Sinks 


Basic  Equation: 


.  .At 
q-  k  A^ 
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Values  of  Thermal  Conductivity  -  "k" 


Material 


Btu : Gm.  Cal. 

hr.'ft.  •  F./ft.       :     sec. 'cm.  '  C./cm, 


Copper  222.0  0.918  (1) 

ru   N^?n    ^  °'145  °-0006  (2) 

(Monofilament) 


Silicone 
Grease 
Dow  Corning 
#  340 


0.363  0.0015  (3) 


(1)  At  18°C.  -  Handbook  of  Chemistry  and  Physics,  Twenty  Fifth  Edition. 

(2)  At  68  F.  -  Elements  of  Material  Science,  by  VanVlack.   Second  Edition 
by  Addison  Wesley,  1964,  p.  420. 

(3)  Dow  Corning  Corporation,  Bulletin  #  04-027,  September  1963. 

Figure  1  shows  the  cross  sectional  view  of  the  heat  sink  copper  plate 
and  the  nylon/grease  separator  between  the  heat  sink  and  the  thermoelectric 
module.   The  total  area  considered  in  these  calculations  is  one  square  foot, 
100  mesh  nylon  material  has  an  open  area  of  approximately  36  percent. 
Therefore  the  nylon  occupies  64  percent  of  the  area  and  the  silicone  grease 
fills  the  remaining  36  percent. 

Assuming  a  ten  degree  F.  temperature  difference  across  the  nylon  and 
grease  separator  the  following  calculations  are  considered: 

thermal      percent  of   temperature 

,         _    ,   conductivity area difference 

heat  transferred  =  — „.  .  . t~~z : : " 

thickness  of  nylon/grease  separator 


Tobler,  Ernst  &  Traber,  Inc.,  Price  List  N-7,  May  1962. 
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COPPER      PLATE 


7 


YzzzzzzzzzzzzzzzMmm 


NYLON 


64% 


GREASE 
*k — 36% A 


—.0208    FT. 


—  .000328     FT. 


Fig.  1.   Cross  section  of  heat  sink  copper  plate 
and  nylon/grease  separator. 
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or: 

=  <kM  A?  W 

^       Ax 

where  Ax   is  determined  as  follows; 

100  mesh  nylon  monofilament  is  approximately  100  microns   in 
diameter.   And, 

1  micron  (10  )  =  1  meter  and   1  meter  =  39.37  inches. 

-4 

_4         qo   )  (39  37) 

Therefore,  100  microns  =  10   meters  =  J ~~ — * — 

Ax   =  0.000328  feet 
ng 

therefore;  q(Nylon)  =  (°' ^OOO^  ^  =  2829  Btu*/Hr'  (10) 

and;      q(Grease)  =  ^63^0^  (10)  .  3984  Btu#/Hr>  (n) 

q(nylon)  +  q(grease)  =  6813  =  -*fc— -  =  ^^^f         (12) 

ng 

therefore;  k   =  0.2235  Btu./hr.*  F.'Ft.   per  foot. 

ng  r 

„  j      o/   i   /      \  -  Ax(nvlon/grease)  m   0.000328  .... 

and;     R(nylon/grease)  =  — *— * — =-J& u   ■  (13) 

Kng 

=  0.00147  ft.2'hr.«  °F./Btu. 

d/                   i         n        Ax(copper  plate)        0.0208  .... 

R(copper  plate)   =  — * — *£ K fc  =  —^ —  (14) 

cu. 
=   0.0000937   ft.2'hr.*°F./Btu. 


Tobler,  Ernst  &  Traber,  Inc.,  Price  List  N-7,  May  1962. 
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R  =   0.00147  +  0.0000937  =  0.0015637  (15) 

and 

C  =  h  =  »  m?CA»  =  639.5  Btu./ft.2'hr.'°F.  (16) 

R       0.0015637 

therefore; 

q  =  C  A  At  =  (639.5)  (0.0273)  (At)  =  17.46  At  Btu./hr.)    (17) 
Equation  (17)  is  shown  graphically  in  Fig.  2. 

Theoretical  Flat  Plate  Heat  Transfer  Calculations 

Empirical  Equation: 


Nu    k  Re     Pr 


Calculation  of  Reynolds  Number: 


L  v 
N   =  -~r         where;  L  =  length  of  plate  in  feet 

v  =  velocity  of  fluid  in  ft./hr. 

2 
V  =  kinematic  viscosity  in  ft.  /hr. 

Given: Y  =  0.0298   ft.2/hr.    @  90°F. 

L  =  0.219  ft.,  length  of  plate  (see  Fig.  3). 

m  =  10  lbs./hr.,  fluid  flow  rate. 

o 
t  =  90  F.  water  inlet  temperature. 

2 
A  =  0.0026  ft.  ,  cross-sectional  area  (see  Fig.  3). 

P   =  62.1  lb  /ft.  ,  water  density  @  90°F. 
m 

k  =  0.354  Btu./hr.  'ft.^F. ,  @  90°F. 


(18) 


For  forced  convection,  laminar  flow. 

** 

Sample  calculations. 
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Fig.  3.   Dimensions  used  for  theoretical  flat  plate 
heat  transfer  calculations. 
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then; 

Volume  of  water  ..   .,  ,  _ 

.,  _  ,,  ..     lbs./hr. 10  ,  1Q* 

flowing       -  v  -  -  -  rrrrrr  ,,0  . :    c  19 ; 

£^  3,  /o^o^\  /iu  /r-  3    (3600)  (62.1) 

ft.  /sec.  (3600)  (lb  /ft. 

IQ 

V  =  0.00004477  ft.  /sec. 
and; 

Velocity  of  ft.3/sec.  (3600)  _  0.00004477  (3600)  ,__. 

water  tlowing   =  v  =  — •— « *— J r4-  =  -1 nAoA  (20) 

-„  ,,  cross-sectional         0.0026 

ft./hr. 

area 

v  =  61.9  ft./hr. 

then; 

.  L^a  (0,219)  (61.9)  =  ( 

\e        V        0.0298        ^  '*  ^ 

and; 

N_  =  5.23 
Pr 

** 
Calculation  of  Heat  Transfer  Coefficient: 

NM  =  0.664  (N_  )*  (N_  )1/3  (22) 

Nu  Re     Pr 

=  0.664  (454.9)^  (5.23)1/3  =  0.664 

=  24.59 
then; 

h  =  CV)  (k>  =  (24.59)  (0.354) 

L  0.219  (23) 

h  =  39.75  Btu./hr.-ft.2'°F. 


** 

Sample  calculations. 
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The  theoretical  values  of  the  heat  transfer  coefficient  are  shown  in 
Table  1.   These  values  are  shown  graphically  in  Fig.  4. 


Theoretical  Heat  Transfer  Data  Plotted  on 
Logarithmic  Scales 


The  data  obtained  from  the  theoretical  analysis  of  flat  plate  heat 
transfer  were  plotted  on  logarithmic  scales.   Graphs  of  equal  water  tempera- 
tures (60  F. ,  75  F. ,  90  F. )  were  obtained  by  plotting  the  heat  transfer 
coefficient  versus  the  mass  flow  rate.   These  data  plotted  as  straight  lines 
on  logarithmic  coordinates  as  shown  in  Fig.  5. 

From  this  data  equations  for  each  water  temperature  line  were  obtained 
as  follows; 


Equation  is  of  the  form:   y  =  bx 


m 


Determination  of  the  slope: 


m  - 


log  y2  -  log  yl 
log  x2  -  log  xl 


=    log  68.83   -   log  39,75 
01  "  log  30  -   log   10 


from  calculations: 


y.  =  39.75  and  x  =  10 


y2  =  68.83 


x2  =  30 


=  1.83778  -  1,59934  = 

m        1.47712  -  i.o        u.*yy* 


(24) 


Determination  of  "y"  Intercept 


y  =   bx 


ra 


(25) 


Therefore; 


b  =  —*•     or     log  b 
m  ° 


=    log  y   -  mlog  x 


And   for  x  =    10 
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h=  (.0383  twa  +  9. 1 325 )  W499>7. 
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WATER  FLOW-"W-LBS./HR. 

h"  versus  mass  flow  rate  "W"  for 
water  temperatures  of  60°F. ,  75°F.  and  90°F.  on  logarithmic  coor- 
dinates. 


Plot  of  heat  transfer  coefficient 

-vOc       tcO 
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log  b  =  log  39.75  -  0.49974  log  10 

log  b  =  1.59934  -  0.49974 

log  b  =  1.09960 
Therefore:    b  =  12.58  ("y"  intercept) 

Equation  of  the  line,  heat  transfer  coefficient  vs.  flow  rate  for  water 
temperature  of  90  F.  is: 

_  ..  QQ  0.4997 
y  =  12.58  x 

Similarly: 

Equation  of  the  line,  heat  transfer  vs.  flow  rate  for  a  water  temperature 

of  60°F.  is: 

.  ..  ,,  0.4997 
y  =  11.43  x 


Determination  of  General  Equation  for  Theoretical 
Heat  Transfer  Coefficient 


From  previous  calculations  it  was  determined  that  the  slope  of  each 
water  temperature  line  was  approximately  equal.   The  temperature  of  the 
water  caused  the  lines  to  move  vertically  on  the  graph. 

The  calculations  also  showed  that  they  intercept  changes  in  equal  seg- 
ments for  each  equal  segment  of  change  in  water  temperature.  Therefore  the 
following  relationship  was  noted: 

For  90  water  temperature  b.  (y  intercept)  is  12.58. 
For  60  water  temperature  b   (y  intercept)  is  11.43. 

and;  Ab  =  b  -  b  =  12.58  -  11.43  =  1.15 

Dividing  Ab  by  30  yields  a  value  of  "b"  for  each  degree  change  in  water 
temperature. 
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Therefore  b.  0p  would  be  11.4683;  b,   0  would  be  11.5066,  etc.   If  the 

above  were  true  from  0  F  to  60  F  then  multiplying  the  water  temperature  by 

.0383  would  yield  the  correct  b  (y  intercept).   However  the  following  was 

noted: 

(60)  (.0383)  =  2.2980   (for  60°F  water  temperature) 

Subtracting  2.2980  from  11.43  (b,^0lJ  =  9.132 

OU  r 

Similarly: 

(90)  (.0383)  =  3.4470   (for  90°F  water  temperature) 
Subtracting  3.4470  from  12.58  (b  op)  =  9.133 
Therefore  the  constant  9.1325  must  be  added  to  each  water  temperature 
multiplied  by  0.0383. 

The  mathematical  relationship  was  determined  to  be; 

b  =  (0.0383)  (water  temperature)  +  9.1325 
The  general  equation  used  to  relate  water  flow,  average  water  tempera- 
ture and  the  heat  transfer  coefficient  was; 

4QQ7 
y  =  (0.0382  t   +  9.1325)  x  '  **'  (26) 

wa 

where; 

y  =  heat  transfer  coefficient  -  Btu./hr.ft.    F. 

o 

t   =  average  water  temperature  -  F. 

x  =  water  flow  rate  -  lbs./hr. 
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Theoretical  Dimensional  Analysis  Calculations 
Assumptions  ( 77   Constant): 


Twic  =  90°F. 


At   =  2°F. 

c 


-rr         _  Twic  _  90  _  ,  c  .        ... 
If ,      =  r- —  =  —  =  45  (constant) 

4    At      2 
c 

W 

7f .     -     -f-   =  0.5  (constant) 
3     w, 
a 

W   =   10  lbs./hr. 
c 

Total  cooling  capacity; 

Q  =  (W  )  (At  )  =  (10)  (2)  =  20  Btu./hr.  (27) 

c     c     c 

Cooling  capacity  per  couple; 

c  20 

%   =  (12)  (3.413)  =  (12)  (3.413)  =  °'488  Wa"S  (28) 

Outlet  water  temperature  on  cooling  side; 

Twoc  =  Twic  -  At  =88  °F.  (29) 

c 

Temperature  of  plate  on  cold  side; 

From  previous  heat  transfer  relationships  it  was  shown  that; 

Qc  =  h  A  Atln   or   h  =  — -  (30) 

In 


Also, 


4997 
y  =  (.0383  t   +  9.1325)  x  '     (Equation  26)  (31) 
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Also  it  is  known  that, 

h  =  y 
Therefore, 
Q 


-=  (.0383  t   +  9.1325)  x  ,4"7  (32) 


A  At,  wa 

In 


where , 


At,         -* 


In     (Twic  -  t   ) 
,  pc 

(Twoc  -  t   ) 
pc 


and, 

x  =  W   ;   t   =  t   +  At,    ;  Twoc  =  Twic  -  At 
c      wa    pc     In  c 

therefore  by  trial  and  error  substitution  for  t   equation  (32) and  be 

pc 

solved  and  the  solution  is; 

o  * 
t   =  70.72  F. 
pc 

Temperature  gradient  across  nylon,  grease  and  copper  plate; 

Q 

At       =  TT7I  =  7^TZ=   1.1455  °F.  (33) 

eg        17.46        17.46 

Temperature  of  cold  junction; 

t  r  =  t   -  At   =  69.5745  °F.*  (34) 

cf    pc     eg 

t   =  5/9  (t  ,  -  32)  =  20.8747  °C 
cc        cf 


IBM  1620  computer  data. 
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Temperature  difference  between  hot  and  cold  junctions; 

1.13  +  .00678(t   )  -  q 

At.   =  n._,   CC ~  =  44.5001  °C.  (35) 

he  .0176 

Temperature  of  hot  junctions; 

t,   ■  t   +  At,   =  65.3748  °C.  (36) 

he    cc     he 

t.  -  =  (9/5)  (t,  )  +  32  =  149.6747  °F.* 
hf  he 

Millivolt  drop  across  each  couple; 

V   =  .526  (At.  )  +  43.4  =  66.8071  MV.  (37) 

ra  he 

Total  power  required  for  module; 


[V  +  (.017)  (I)  (t   )]  [(12)  (I)] 


• 


P   =  — 1„nnCC =  12.9834  watts   (38) 

w  1000 

P  .  -  (P  )  (3.413)  =  44.3124  Btu./hr.* 
wo     w 

Total  heat  to  be  rejected  at  hot  junctions; 

Qu     =   Q     +PU=20+  44.3124       64.3124  Btu./hr.*  (39) 

h  xc  wb 

Temperature  gradient  across  nylon,  grease  and  copper  plate; 

At   =  t—tt  =  3.6834  °F.  (40) 

eg    17.46 

Temperature  of  plate  on  hot  side; 

t  .  =  t.  .  -  At   =  145.9913  °F.*  (41) 

ph    hf     eg 

Water  flow  rate  through  hot  side; 

Wc    10 
Wh  =  ^  =  -^  =  20  lbs./hr.  (42) 
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Temperature  rise  of  water  through  hot  side; 

h         o  * 
Ath  =  —  =  3.2156  °F.  (43) 

h 

Inlet  water  temperature  at  hot  side; 

From  previous  heat  transfer  relationships, 

Qh  =  h  A  At     or   h  =  j-jj-  (44) 

In 


and, 


4997 
y  =  (.0383  t   +  9.1325)  x       (Equation  26)       (45) 
wa 


and  it  is  known  that, 

h  =  y 
where , 

At 


In     (Twih  -  t  ,  ) 
,  ph 

(Twoh  -  t  .  ) 
ph 


and, 


x  =  W,    ;    t   =  t  .  -  At. 
h        wa    ph     In 


therefore, 

Q 


h  4997 

-=  (.0383  t   +  9.1325)  x  '  **'  (46) 


A  At,  wa 

In 


and  by  trial  and  error  substitution  for  Twih  equation  (46) can 
be  solved  and  the  solution  is, 

Twih  =  104.7413  °F.* 
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Outlet  water  temperature  at  hot  side; 

Twoh  =  Twih  +  At  =  107.9769  °F.  (47) 

n 

77  =  — £  =  — 2S =  4513*  (48) 

wb 

^   .  Twih^-  Twic  ,  104,7413  -  90,  ,  7>37Q7*  ^ 

c 

W 

T->   =  7T  =  ™  =  '3  (constant)  (50) 

3   Wh   20 


77/  =  77 —  =  —  =  45  (constant)  (51) 

c 


Assumptions  ~ff     Constant: 


Twic  =  90°  F. 


At   =  2°  F. 
c 


-rr  TwJC     90     .  _   . 

II  .      =  T~ —  =  ~  =  45  (constant) 


Twih  =  110°  F. 


_.   _  Twih  -  Twic  _  110  -  90  _  ..  ,  ... 

772  =  — =  r 10,  ^constant) 

c 


W   =10  lbs./hr. 
c 

A   =  .0273  ft.2 


Calculations: 


Total  cooling  capacity; 

Q  =  (W  )  (At  )  =  (10)  (2)  =  20  Btu./hr.  (52) 

c     c     c 
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Cooling  capacity  per  couple; 


=  0.488  watts  (53) 


Nc   (12)  (3.413) 

Outlet  water  temperature  on  cooling  side; 

Twoc  =  Twic  -  At  =  90  -  2  =  88°  F.  (54) 

c 

Temperature  of  plate  on  cold  side; 

From  previous  heat  transfer  relationships  it  was  shown  that, 

Qc  =  h  A  At     or   h  =  —-  (55) 

In 


Also, 


4997 
y  =  (.0383  t   +  9.132)  x  *     (Equation  26)         (56) 

J  wa  ^ 


Also  it  is  known  that, 

hSy 
Therefore, 
Q. 


=  (.0383  t   +  9.132)  x  ,4"7  (57) 


A  At,  wa 

In 


where, 


At.         -*- 


In     (Twic  -  t   ) 

,  PC 

1    (Twoc  -  t   ) 
pc 


and, 

x  =  W    ;    t   =  t   +  At,     ;   Twoc  =  Twic  -  At 
c   '    wa    pc     In  < 
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therefore  by  trial  and  error  substitution  for  t   equation  (  57)  can 

PC 

be  solved  and  the  solution  is; 

t   =  70.72  °F.* 
pc 

Temperature  gradient  across  nylon,  grease  and  copper  plate; 

Q 

At   =  T7TZ  =  7TZZ  =  l«1455  °F-  (58) 

eg    17. 46    17.46 

Temperature  of  cold  junctions; 

t  .  =  t   -  At   =  69.5745  °F.*  (59) 

cf    pc     eg 

t   =  5/9  (t  ,  -  32)  =  20.8747  °C.  (60) 

cc        cf 

Temperature  difference  between  hot  and  cold  junctions; 

1.13  +  .00678  (t   )  -  q  ^ 

At   =  cc c  m   44#5001  °c#  (61) 

he  .0176 

Temperature  of  hot  junction; 

t.      ■  t   +  At.   =  65.3748  °C.  (62) 

he    cc     he 

t,  .  =  (9/5)  (t.  )  +  32  =  149.6747  °F.* 
hf  he 

Millivolt  drop  across  each  couple; 

V   =  .526  (At,  )  +  43.4  =  66.8071  MV.  (63) 

m  he 

Total  power  required  for  module; 


[V  +  (.017)  (I)  (t   )]  [(12)  (I)] 


• 


P  =  ~ ,nnn  CC =  12.9834  watts    (64) 

w  1000 


* 
IBM  1620  computer  data. 
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P  ,  =  (P  )  (3.413)  =  44.3124  Btu./hr.*  (65) 

wb     w 

Total  heat  to  be  rejected  at  hot  junctions; 

Qu  =  Q  +  p  u  ~  64.3124  Btu./hr.*  (66) 

h     c    wb 

Temperature  gradient  across  nylon,  grease  and  copper  plate; 

Q, 

At   =  777 t  =  3.6834  °F.  (67) 

eg   17.46 

Temperature  of  plate  on  hot  side; 

t  ,  =  t,  .  -  At   =  145.9913  °F.*  (68) 

ph    hf     eg 

Water  flow  through  hot  side; 

From  previous  heat  transfer  relationships; 

Qh 
Q  =  h  A  At,    or   h  =  ;   "  (69) 

xh         In  A  At, 

In 


and, 


4997 

y  =  (.0383  t   +  9.132)  x  *  **'  (70) 

J  wa 


and  it  is  known  that, 

h  =  y 
where, 

At, 


At 


'h 


T   ;   x  =  W,    ;   t   =  t  .  -  At. 
)  h       wa    ph     In 


In     (Twih  -  t  .  )   '  h       wa    ph 

,  ph 

in(Twoh  -  t  .  ) 
ph 


and; 
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^h  4997 

7—^ —  =  (.0383  t   +  9.132)  x  *   *'  (71) 

A  At,  wa 

In 

therefore  by  trial  and  error  substitution  for  w,  equation  (71 ) 
can  be  solved  and  the  solution  is; 

Wh  =  25.35  lbs./hr. 

Temperature  rise  of  water  through  hot  side; 

Qh        _  * 
Ath  =  ~  =  2.5370   F.  (72) 

h 

Outlet  water  temperature  at  hot  side; 

Twoh  =  Twih  +  At,  =  112.5370  °F.  (73) 

n 

Therefore; 

Q 
T    =  =-*  =  .4513   *(variable)  (74) 

wb 


^2  =  TWlhAt  TWi°  =  10*°   *<constant>  (75) 


W 
7",  *  —■  ■  0.3945   *(variable)  (76) 

J   Wh 


7/*4  =  —^  =  45   *(constant)  (77) 

c 
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Analysis  of  Theoretical  Data 


The  theoretical  data  acquired  from  the  IBM  1620  computer  were  pro- 
grammed for  multiple-regression  analysis  by  an  IBM  1410  computer.   Plates 
X  and  II  graphically  show  the  theoretical  data*  Appendix  C  shows  the  data 
furnished  the  IBM  1410  computer. 

The  following  equations  were  used  for  the  statistical  analysis: 


=  A  +  B  T    +  C  7T  +  D  J. 
2      3      4 


=  A  +  B  7T-  +  C  J-  +   D  7T,    +  E  7T0  7"0  T. 
2      3      4      2  3  4 

=  A  +   B72+   C  7f     +   D7"4+   E  7T2  T^  7Th  +   FF4' 


1. 

If 

2. 

ri 

3. 

Wl 

4. 

Wl 

5. 

7Tl 

6. 

7Tl 

7. 

wl 

=  a  +  b  w2  +  c  j3  +  d  ^2  r3  irk 


=  A+B^2+C73  +  D74' 


=  A  +   B  77*0  +   C  7T_  +  D  /7".    +   E  7"/ 
2  3  4  4 


2  3  2  3  4  4 


The  data  also  were  converted  to  natural  logarithmic  values.   These 
values  were  also  processed  by  the  IBM  1410  computer  for  the  best  fit  analy- 
sis of  an  equation  of  the  form; 


log  77^  =  A .+  BlogT2  +  Clog7*3  +  Dlog7"4 


(78) 


From  equation  (78)  an  equation  of  exponential  form  can  be  derived  as 
follows; 


77*,  =  eA  (/7\)B  (7T  )C  {7T,)^ 
1         2      3      4 


(79) 


From  the  multiple  regression  analysis  by  the  IBM  1410  computer  the 
following  polynomial  equation  was  determined  to  provide  adequate  precision 


EXPLANATION  OF  PLATE  I 


Plot  of  the  computer  data  of  the  pi  term  Q  /P  .  on  the  pi  term 

CT  .,  -  T  .  )/At  while  holding  the  pi  term  W  /W,  constant  and 
N  wih    wic   c  °     r        ch 


the  pi  term  T  .  /At  constant  at  four  values, 
wic   c 
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EXPLANATION  OF  PLATE  II 


Plot  of  the  computer  data  of  the  pi  term  Q  /P   on  the  pi  term 

W  /W,  ,  while  holding  the  pi  term  (T  . ,  -  T  .  VAt  constant  and 
c  h  o     r      \  W1h    W1C   c 

the  pi  term  T  .  /At  constant  at  four  values. 

r  W1C    c 
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for  engineering  design; 

7fl   =   0.55162   -  0.02902  7T     -  0.56858  7^  +   0.00752/7^  (80) 

(R2  =   0.996    ,    4>         =   0.00945) 
yx 

The  exponential  equation  obtained  from  the  computer  was  as  follows; 

Tx   =  0.091  (7T2  ♦  ID)"0'4333  C^)-0-3945  (^4)°-6454      (81) 

(R2  =  0.835  ,  ^         =  0.13641) 
yx 

Since  77~_  had  negative  as  well  as  positive  values,  the  data  for  the 
logarithmic  form  had  to  be  converted  to  positive  values  by  adding  a  factor 
of  10  to  each  value  of  7T   ,      This  explains  the  reason  for  (  7"^  +  10) 
appearing  in  the  final  equation. 

The  1410  multiple  regression  computer  program  is  not  included  in  the 
thesis  since  it  is  a  common  item  furnished  by  the  computing  center. 


EXPERIMENTAL  EQUIPMENT  AND  INSTRUMENTATION 

The  equipment  used  in  this  study  consisted  of  the  following: 

1.  A  twelve-couple  thermoelectric  module  of  the  type  now  commercially 
available. 

2.  Two  flat-plate  design  heat  sinks  equipped  for  liquid  flow. 

3.  A  2-foot  square  wooden  box  for  housing  the  device  and  to  provide 
for  insulation  of  the  device. 

4.  Two  precisely  controlled  water  baths,  with  stirrers,  thermo- 
regulators  and  overflow  provision.  The  overflow  was  used  to 
maintain  a  constant  head  on  the  device. 

5.  A  vacuum  tank  for  removing  air  from  the  distilled  water. 
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6.  A  return  tank  and  a  supply  tank  with  suitable  pump. 

7.  A  d.c.  power  supply  consisting  of  two  six  volt  automobile  type 
batteries  connected  in  parallel  and  two  rheostats  to  control  the 
current. 

8.  A  thermocouple  switch. 

9.  An  electronic  ice  bath. 

10.  A  precision  potentiometer. 

11.  An  ammeter  and  voltmeter. 

12.  A  precision  balance  to  batch  weigh  the  liquid. 

13.  Needle  valves  to  control  flow  rate. 

14.  Liquid  filters  and  associated  apparatus. 

15.  Calibrated  thermocouples  for  determining  the  water  temperatures 
and  degrees  of  cooling  and  heating  of  the  water  by  the  device. 

16.  Thermocouples  soldered  to  junctions  to  determine  hot  and  cold 
junction  temperatures. 

The  liquid  heat  sinks  were  made  of  copper.   The  heat  sinks  were  approx- 
imately 1.5  inches  wide  by  2.25  inches  long  so  as  to  fit  the  heat  transfer 
area  of  the  thermoelectric  module. 

A  mixing  tube  was  attached  to  the  outlet  of  each  sink.  The  purpose  of 
the  mixing  tube  was  to  assure  adequate  mixing  of  the  liquid  before  the  out- 
let thermocouple  sensed  the  liquid  temperature.  The  mixing  tube  was  made 
from  one-half  inch  diameter  copper  tubing.  Inside  the  tube  and  at  one  end 
a  small  copper  disc  was  suspended  by  means  of  thin  spiders  in  the  center  of 
the  tube.  Liquid  flowing  through  the  tube  had  to  pass  over  and  around  this 
disc  thus  causing  turbulence  and  subsequent  mixing  of  the  liquid.  A  cross- 
section  of  the  device  is  shown  in  Plate  III.   Plate  IV  is  a  photograph  of 


EXPLANATION  OF  PLATE  III 

Cross-section  of  thermoelectric  heating/cooling  device  showing 
arrangement  of  thermoelectric  module,  heat  sinks  inlets  and  out- 
lets, thermocouples  and  associated  equipment. 
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EXPLANATION  OF  PLATE  IV 

View  of  thermoelectric  heating/cooling  device  mounted  inside 
plywood  container  and  ready  for  addition  of  insulation 
material. 


PLATE   IV 
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the  device  mounted  inside  the  plywood  container  and  ready  for  addition  of 
the  insulating  material. 

Short  lengths  of  plastic  tubing  were  used  to  connect  the  inlets  and 
outlets  to  3/8  inch  diameter  copper  tubing  which  connected  to  the  water  bath 
outlets. 

The  block  diagram  of  the  laboratory  equipment  is  shown  in  Plate  V.   The 
water  baths  were  equipped  with  overflow  outlets  which  provided  the  device 
with  a  constant  head  of  water.   By  this  system  of  overflows  precise  rates 
of  flow  were  possible  to  attain  and  hold  constant.   Each  water  bath  con- 
tained an  electric  stirrer  and  a  precision  mercury  actuated  thermoregulator. 
The  mercury  actuater  was  electrically  connected  to  a  solid-state,  low  cur- 
rent actuated,  5-ampere  capacity  relay.   This  relay  controlled  the  power  to 
the  300-watt  low  time  lag  water  heaters  in  each  water  bath.   Each  water  bath 
outlet  was  provided  with  a  fiberglass  water  filter.   Plate  VI  is  a  photograph 
of  the  experimental  system  as  it  appeared  in  the  laboratory. 

A  twenty  gallon  plastic  supply  tank  was  provided  to  furnish  the  dis- 
tilled water  for  both  water  baths.   Water  could  be  circulated  through  a 
heat  exchanger  before  going  into  the  water  bath  is  a  lower  water  bath  temper- 
ature was  required.   Water  entering  the  water  baths  was  maintained  several 
degrees  cooler  than  the  desired  water  bath  temperature.   The  water  bath 
heaters  then  heated  the  water  to  the  desired  temperature  and  turned  off  as 
determined  by  the  thermoregulators.   Cooling  water  was  obtained  from  a 
standard  electric  water  fountain  type  cooler. 

A  vacuum  tank  was  constructed  from  a  liquified  gas  cylinder.  A  vacuum 
pump  was  used  to  obtain  a  high  vacuum  (about  28  inches  Hg. )  on  a  portion  of 
the  distilled  water  which  was  in  the  tank.   The  vacuum  pump  was  allowed  to 


EXPLANATION  OF  PLATE  V 


Block  diagram  of  the  Laboratory  equipment. 
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EXPLANATION  OF  PLATE  VI 


View  of  the  experimental  equipment  in  the  laboratory, 


PLATE  VI 
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operate  for  4  to  6  hours  to  "boil  off"  the  air  from  the  water.   This  deaer- 
ated  water  was  then  removed  from  the  vacuum  tank  by  means  of  a  water  syphon 
system  and  then  pumped  into  the  main  supply  tank  for  release  into  the  water 
baths.   It  was  found  that  deaerated  water  helped  to  maintain  a  more  constant 
flow  through  the  system.   Air  in  the  water  can  accumulate  at  various  loca- 
tions in  the  system  particularly  in  the  needle  valves  and  precision  flow- 
raters  thus  causing  variations  in  the  water  flow  rate. 

A  return  tank  provided  a  method  of  collecting  the  liquid  from  the  over- 
flow tubes  and  the  outlets  from  the  thermoelectric  device.   A  pump  connected 
to  the  return  tank  allowed  the  liquid  to  be  either  transferred  to  the  vacuum 
tank  for  deaerating  or  to  the  main  supply  tank  for  return  to  the  water 
baths.   A  high  quality  water  filter  was  installed  after  the  water  pump  to 
further  assure  clean  water  for  the  system. 

A  500-gram  precision  scale  provide  for  accurate  batch  weighing  of  300 
gram  samples  of  the  liquid  flowing  through  each  circuit.   Electric  solonoid 
valves  provided  for  selection  of  the  liquid  from  each  flow  circuit  for  batch 
weighing.   A  sample  was  taken  each  time  temperature  and  power  data  were 
noted.   Any  change  in  flow  rate  was  thus  quickly  ascertained  and  corrected. 

Precision  flowrators  were  used  to  determine  the  approximate  rate  of 
flow  through  each  circuit.  However  final  determination  of  the  flow  rate 
was  made  with  the  precision  scale. 

A  diagram  of  the  temperature  measurement  instrumentation  is  shown  in 
Plate  VII.   Temperatures  were  measured  using  calibrated  copper-constantan 
thermocouples  connected  to  a  high  quality  thermocouple  switch.   The  output 
terminals  on  this  switch  were  connected  to  the  input  terminals  of  an  elec- 
tronic ice  bath.   The  electronic  ice  bath  provided  a  highly  stable  reference 


EXPLANATION  OF  PLATE  VII 


Diagram  of  temperature  measurement  instrumentation. 
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junction  for  the  measurement  system  without  the  need  for  inaccurate  water 
type  ice  baths.   The  output  from  the  electronic  ice  bath  provided  a  precise 
millivolt  signal  that  could  be  accurately  measured  by  a  precision  potenti- 
ometer.  The  supply  voltage  to  the  electronic  ice  bath  was  precisely 
regulated  to  further  assure  accurate  reference  junction  output. 

The  difference  between  the  inlet  and  outlet  water  temperatures  was 
determined  by  connecting  the  inlet  and  outlet  thermocouples  so  as  to  buck 
each  other.   This  connection  was  provided  at  the  thermocouple  switch  by 
connecting  the  constantan  terminals  together  with  a  short  length  of  constan- 
tan  wire.   The  output  from  these  thermocouples  was  then  connected  to  a 
double  pole-double  throw  knife  switch  and  thence  to  the  input  terminals 
on  the  potentiometer.   When  it  was  desired  to  measure  the  temperature  dif- 
ference the  thermocouple  switch  was  placed  in  position  zero  (all  contacts 
open).   The  knife  switch  could  then  be  positioned  so  as  to  measure  the 
temperature  difference.   The  electronic  ice  bath  is  not  connected  into  the 
measurement  circuit  during  measurement  of  the  temperature  difference, 
because  of  the  open  contacts  in  the  thermocouple  switch. 

Plate  VIII  is  a  diagram  of  the  direct  current  power  supply  for  the 
thermoelectric  module  and  the  electrical  measurement  instrumentation.   The 
direct  current  power  for  the  thermoelectric  module  was  obtained  from  two 
6-volt  automotive  type  lead  acid  batteries  connected  in  parallel.   Two 
rheostats  were  connected  in  parallel  and  inserted  in  one  of  the  power  lines 
to  the  module.   One  rheostat  was  used  to  control  the  largest  portion  of  the 
voltage  drop  to  the  module  from  the  batteries.   The  second  rheostat  was  used 
to  drop  a  small  portion  of  the  voltage  and  acted  as  a  fine  adjustment  for 
the  current  to  the  module.   By  this  means  the  current  to  the  module  could 


EXPLANATION  OF  PLATE  VIII 

Diagram  of  direct  current  power  supply  for  thermoelectric  module 
and  electrical  instrumentation. 
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be  precisely  controlled. 

A  direct  current  ammeter  and  voltmeter  were  used  to  measure  the  power 
to  the  module.   The  ammeter  had  a  full  range  of  15  amperes  and  the  voltmeter 
a  full  scale  range  of  one  volt.   The  voltmeter  could  be  read  to  .005  volts. 

EXPERIMENTAL  DATA 

Analysis  of  Laboratory  Data 

The  laboratory  data  were  programmed  for  multiple-regression  analysis 
by  an  IBM  1410  computer.   Plates  IX  and  X  graphically  show  the  laboratory 
data  along  with  the  theoretical  data.  Appendix  C  shows  the  data  furnished 
the  IBM  1410  computer. 

The  following  equations  were  used  for  the  statistical  analysis: 

1.  77"   =  a  +  b^  +  cT3 

2.  7Fl  =  a  +  B7T2  +  c7T3  +  d^27"3 

3.  7Tl  =  A  +  B  7T2  +  C  7T3  +  D  7T2  T3  +  E  7T 3 

4.  7Tl  =  A  +  B  7T2  +  c  7"3  +  D  7] 32 

The  data  were  also  converted  to  natural  logarithmic  values.   These 
values  were  also  processed  by  the  IBM  1410  computer  for  the  best  fit  analy- 
sis of  an  equation  of  the  form; 

log  7Tl  =   A  +  Blog  7T2  +   Clog  7T3  (82) 

from  equation  (82)  an  equation  of  exponential  form  can  be  derived 
as  follows; 

Tx  =  *A  <>72B)  (7T3C)  (83) 


EXPLANATION  OF  PLATE  IX 


Plot  of  theoretical  and  laboratory  data  of  the  pi  term  Q  /P   on 

the  pi  term(T  ,. -T  .  )/At  while  holding  the  pi  term  W  /W,  constant 
r  wih  wic   c  c  c   h 

and  the  pi  terra  T  .  /At  constant  at  selected  values, 
wic   c 
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EXPLANATION  OF  PLATE  X 


Plot  of  theoretical  and  laboratory  data  of  the  pi  term  Q  /P  ,  on 

J  r        x   wb 

the  pi  term  W  /W,  while  holding  the  pi  term(T  . ,  -T  .  )/At  constant 
r        c  h  r      \  W1h  W1C   c 

and  the  pi  term  T  .  /At  constant  at  selected  values. 
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From  the  multiple  regression  analysis  by  the  IBM  1410  computer  the 
following  polynomial  equation  was  determined  to  provide  adequate  precision 
for  engineering  design; 

7fl   =  1.18418  -  0.02457  W 2   -   0.05885  ^~3  (84) 

(R2  =  0.859  ,  >    =  0.02693) 
yx 

The  exponential  equation  obtained  from  the  computer  was  as  follows; 

T1   =  2.094!  (  7T2)-°-37515  (7r3,-°-°"66  (85) 

(R2  =  0.874  ,  A>        =  0.03015) 
y«x 


Comparison  of  Experimental  and  Theoretical  Heat 
Sink  Heat  Transfer  Coefficients 


From  the  comparison  of  laboratory  data  with  theoretical  data  a  large 
difference  in  coefficient  of  performance  was  noted.   It  was  concluded  that 
part  of  this  difference  could  be  accounted  for  by  a  difference  in  the  heat 
transfer  coefficient  of  the  laboratory  device  as  compared  to  a  theoretical 
flat  plate. 

A  comparison  can  be  derived  by  using  the  junction  temperatures  as  a 
base  for  calculating  the  log  mean  temperature  difference  between  the  junc- 
tions and  the  water  being  cooled  or  heated. 

The  equation  used  was  as  follows: 

h  =  -a- 


AAt. 
In 


where 


h  =  heat  transfer  coefficient  -  Btu/hr  ft   F 
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At 


2 
A  =  area  of  heat  transfer  surface  -  ft 

q  =  rate  of  heat  transfer  -  Btu/hr 

At 
In      t  .  -  t. 

ln  t   -  t 

WO      J 

and: 

t  .  =  inlet  water  temperature  -  F 

o 

t   =  outlet  water  temperature  -  F 

wo  r 

t.   =  junction  temperature  -  F 

The  calculated  heat  transfer  coefficients  are  shown  in  Tables  2  and  3. 
Also  shown  are  the  ratio  between  the  experimental  values  and  the  theoretical 
values.   From  this  data  it  was  concluded  that  the  heat  transfer  coefficient 
of  the  experimental  device  was  greater  than  the  coefficient  for  a  theoretical 
flat  plate  under  the  assumed  conditions.   The  ratio  of  the  experimental  to 
the  theoretical  heat  transfer  coefficient  related  to  the  cold  and  hot  junc- 
tions was  about  3.32:1  and  2.46:1,  respectively. 

DISCUSSION 

The  results  of  this  research  indicate  that  the  methods  and  equipment 
used  to  conduct  the  laboratory  tests  were  adequate  to  determine  dimension- 
less  term  relationships  for  engineering  design  requirements.   A  device  hav- 
ing larger  capacity  for  cooling  could  improve  on  the  accuracy  of  laboratory 
measurements  by  providing  a  greater  temperature  drop  as  the  liquid  passes 
through  the  heat  exchanger.   Also  the  flow  rates  could  then  be  increased  by 
providing  adequate  systems  for  handling  the  flow  of  liquid  and  for  liquid 
bath  heating  and/or  cooling. 
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Multiple  regression  was  used  to  determine  the  relationships  between  the 
pi  terms  involved  in  the  analysis.   Polynomial  equations  of  the  form 
77     =  A  +  B  77     +  C  71     and  exponential  equations  of  the  form  77     =  C  7f      77 
were  developed  from  computer  programs.   The  equation  showing  adequate  fit  to 
the  data  and  having  the  least  number  of  terms  was  chosen  from  the  group  of 
equations  solved  by  the  computer  programs. 

As  shown  in  Plate  IX  the  relationship  of  77".  ( Qc / Pwb )  to  7" 
(Twih-Twic/Atc)  was  almost  linear  with  7T  (Wc/Wh)  held  constant.   Plate  X 
shows  a  non-linear  relationship  of  7T   (Qc/Pwb)  to  77   (Wc/Wh)  with 
77 « (Twih-Twic/Atc)  held  constant.   The  above  relationships  were  approxi- 
mately true  for  both  theoretical  and  experimental  data. 

The  mercury  thermoregulators  used  to  maintain  constant  water  bath 
temperatures  were  very  sensitive.   Difficulty  was  experienced  when  mounting 
these  regulators  directly  on  the  sides  of  the  water  baths.   The  slightest 
vibration  from  the  stirrer  motors  could  cause  the  regulator  relay  to 
chatter.   The  final  solution  to  this  problem  was  to  mount  the  regulators 
on  a  rod  which  was  isolated  from  the  water  baths. 

CONCLUSIONS 

The  following  conclusions  were  drawn  from  the  analysis  and  data 
presented : 

1.   The  method  of  dimensional  analysis  can  be  used  to  determine  the 
operating  characteristics  of  a  thermoelectric  heating/cooling  device.   The 
device  can  be  either  operated  with  liquid  or  air  flowing  through  the  heat 
sinks  or  liquid  through  one  sink  and  air  through  the  other.   For  example  the 
device  could  be  used  to  cool  and  heat  water  or  cool  water  and  heat  air. 
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2.  The  slopes  of  the  graphs  for  the  experimental  and  theoretical  data 
were  not  exactly  the  same.   It  is  believed  that  experimental  error  can 

account  for  part  of  this  inaccuracy.   The  experimental  data  showed  that  for 

■77*  Twi  h  —  Twi  c  ~tt     Wc  /  /* 

//„( 7"7 )  constant,  //  „  (~)  decreased  in  value  as  '/  .  (cooling  CO. P.) 

2  Ate  3  Wh  1        ° 

increased  in  value  and  appeared  to  be  a  curvilinear  function.   Also  when 

''^(777)  was  held  constant,  f/n( 77 )  decreased  as  7T .    (cooling  CO. P.) 

3  Wh  2     Ate  1 

increased  and  the  relationship  was  essentially  linear.  Both  relationships 
were  in  general  agreement  with  the  predicted  functions  as  derived  from  the 
theoretical  analysis. 

3.  Another  factor  which  could  have  caused  a  deviation  between  experi- 
mental and  theoretical  results  could  be  a  difference  in  the  heat  transfer 
coefficient  of  the  experimental  device.   The  theoretical  calculations  were 
based  on  an  emperical  equation  for  a  flat  plate  under  laminar  flow  condi- 
tions.  Although  the  experimental  device  approached  a  flat  plate  configura- 
tion, the  liquids  entered  and  left  the  heat  transfer  system  at  right  angles. 
This  could  have  caused  entrance  and  exit  effects  that  the  theoretical 
analysis  did  not  include. 

4.  The  experimental  data  did  not  show  any  measurable  difference  due  to 
different  water  inlet  temperatures  at  the  cooling  side  of  the  device.   From 
the  theoretical  data  an  error  of  less  than  0.1  degree  Fahrenheit  in  the 
measurement  of  the  temperature  drop  of  the  water  being  cooled  could  have 
nullified  the  theoretical  difference  due  to  the  inlet  water  temperature. 

Since  the  inlet  water  temperature  appeared  to  have  no  significant  effect, 

Twic 
the  fourth  Pi  term  (— - — )  was  eliminated  from  the  statistical  analysis  of 

the  experimental  data. 

5.  An  empirical  equation  for  the  type  of  heat  sinks  used  in  this 
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research  was  not  available.   Therefore  close  agreement  between  actual  values 
of  experimental  and  theoretical  data  could  not  be  expected. 

SUGGESTIONS  FOR  FURTHER  RESEARCH 

Laboratory  experiments  of  the  type  described  here,  require  highly 
accurate  instruments.   Instruments  for  measuring  the  electrical  power  deliv- 
ered to  the  device  should  be  as  accurate  as  possible.   Thermocouples  for 
determining  the  temperature  drop  of  the  liquid  being  cooled  should  be 
calibrated. 

It  was  possible  with  the  instruments  available  and  calibrated  thermo- 
couples to  measure  a  temperature  change  of  2  F  in  the  liquid  being  cooled. 
In  future  research  a  larger  device  should  be  assembled  so  as  to  provide  a 
greater  heat  transfer  capacity.   This  greater  heat  transfer  capacity  would 
result  in  a  greater  temperature  drop  of  the  liquid  being  cooled.   The  result 
should  be  a  greater  degree  of  precision  in  the  measurement  of  the  temperature 
drop  of  the  liquid. 

Other  investigators  who  may  perform  similar  research  should  carefully 
consider  the  precise  control  of  liquid  flow,  and  assure  themselves  of  reli- 
able equipment.   A  closed  system  for  the  liquid  circulation  with  suitable 
filters  should  be  considered  for  future  research.   However  the  open  type 
system  and  filters  used  in  these  tests  proved  to  be  satisfactory. 

Future  research  should  be  carried  out  with  other  temperature  drops  of 
the  liquid  being  cooled.   The  tests  reported  were  performed  for  the  liquid 
being  cooled  2  F.   Additional  research  would  provide  a  family  of  equations 
related  to  various  temperature  drops.   Such  equations  could  then  be  applied 
to  engineering  design  for  a  particular  device. 
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Also  these  tests  were  performed  with  a  current  flow  of  15  amperes. 
Future  tests  should  be  conducted  on  a  device  with  other  currents,  i.e.  20, 
25,  and  30  amperes,  or  up  to  the  limit  of  current  for  the  particular  thermo- 
electric module  being  used. 

Future  research  should  encompass  the  preliminary  testing  of  a  particu- 
lar device  to  determine  the  heat  transfer  coefficient  of  the  heat  sinks,  if 
a  satisfactory  theoretical  heat  transfer  coefficient  is  not  available. 
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APPENDIX  A 

Determination  of  Conversion  Factor  to  Convert 
Thermal  Conductivity  from  Metric  Units  to 
British  Thermal  Units 
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Determination  of  conversion  factor  to  convert  from  gram  calories  per 
square  centimeter  per  degree  Centigrade  per  second/centimeter  to  British 

Thermal  Units  per  square  foot  per  degree  Fahrenheit  per  hour/foot. 

-3 

Gram  calories  x  3.968  x  10  "  =  British  Thermal  Units 

Degrees  Centigrade  x  1.8  =  degrees  Fahrenheit 

-3 

Square  centimeters  x  1.076  x  10   =  square  feet 

Centimeters 

30.48    =  feet 

Seconds 

l6bT~  =  hours 

and; 

gm.  cal. 


k  = 


2  °^  / 
sec.  cm.    C./cm. 


Substitution  of  these  conversion  units  in  the  above  equation  results  in  the 
following  relationship; 

k  _  (gm.  cal.)(3.9681)(10~3) 


(OC.)(1.8)(cm.2)(1.076)(10"3)(^^:)/^a7^ 

3600   30.48 


k  =    (3.968)(3600) 

(1.8)(1.076)(30.48) 

Btu. 


k  =  241.97 


hr.  ft.2  °F./ft. 


Therefore: 

B|u^ m      gm^cj^ (241>97) 

hr.  ft.    F./ft.     sec.  cm.    C./cm. 
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APPENDIX  B 


Sintered  Frigistor  Performance  Data 
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APPENDIX  C 

Data  Programmed  for  Multiple-regression 
Analysis  by  IBM  1410  Computer 
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Table  3.      Condensed   theoretical  data  and   program  data   for 

multiple   regression  analysis  by    IBM  1410   computer, 


THEORETICAL    POLYNOMIAL 


.2152 

10.0 

.2646 

10.0 

.3164 

10.0 

.3707 

10.0 

.4278 

10.0 

.4877 

10.0 

.5508 

10.0 

.2554 

10.0 

.3051 

10.0 

.3571 

10.0 

.4116 

10.0 

.4688 

10.0 

.5288 

10.0 

.5919 

10. c 

.2946 

10.0 

.3445 

10.0 

.3967 

10.0 

.4513 

10.0 

.5086 

10.0 

.5686 

10.0 

.6316 

10.0 

.3328 

10.0 

.3829 

10.0 

.4352 

10.0 

.4899 

10.0 

.5472 

10.0 

.6072 

10.0 

.6701 

10.0 

.6329 

35.0 

.5797 

35.0 

.5036 

35.0 

.4113 

35.0 

.3141 

35.0 

.220 

35.0 

.1355 

3^.0 

.6780 

40.0 

.6009 

40.0 

.5079 

40.0 

.4045 

40. 0 

.3012 

40.0 

.2043 

40.0 

.12  01 

40.0 

.7071 

45.0 

.6107 

45.0 

.5056 

45.0 

.3945 

45.0 

.2861 

45.0 

.1876 

45.0 

.1048 

45.0 

.7240 

50.0 

.6143 

50.0 

.4975 

50.0 

.3800 

50.0 

.2688 

50.0 

.1703 

50.0 

.0898 

5  0.0 

221. 

.5150 

1225 

.0 

202, 

,8950 

1225 

.0 

176 

.260 

1225 

.0 

143, 

.9550 

1225 

.0 

109 

.9350 

1225 

.0 

77, 

.0 

1225 

.0 

47, 

.4250 

1225 

.0 

271 

.20 

1600 

.0 

240 

.360 

1600 

.0 

203 

.160 

1600 

.0 

161 

,80 

1600 

.0 

120 

.480 

1600 

.0 

81 

.720 

1600 

.0 

48, 

.040 

1600 

.0 

318, 

,1950 

2025 

.0 

274 

.8150 

2025 

.0 

227 

.520 

2025 

.0 

177 

.5250 

2025 

.0 

128, 

,7450 

2025 

,0 

84 

.3750 

2025 

,0 

47, 

,160 

2025 

,0 

362, 

.0 

2500 

,0 

307, 

,650 

2500 

.0 

248, 

,750 

2500 

,0 

190, 

.0 

2500, 

.0 

134, 

,40 

2500, 

,0 

85, 

,150 

2500, 

.0 

44, 

.90 

2500, 

,0 

Table   3    (cont. ). 


85 


THEORETICAL  PCLYN 

CM  I  A  L 

** 

r, 

Ts 

7% 

7T4 

T2T377i- 

.2152 

14.  1042 

.50 

35.0 

246.8235 

1225.0 

.2646 

12.3314 

.50 

35.0 

215.7995 

1225.0 

.3164 

10.1093 

.50 

35.0 

176.912a 

1225.0 

.3707 

7.5992 

.50 

35.0 

132.9860 

1225.0 

.4278 

4.8803 

.50 

35.0 

85.4053 

1225.0 

.4877 

2=*527 

.50 

35.0 

35.9223 

1225.0 

.5  508 

-0.8905 

.50 

35.0 

-  15.5838 

1225.0 

.6172 

-  3.9311 

.50 

35.0 

-68.7943 

1225.0 

.6873 

-7.0192 

.50 

35.0 

-122.8360 

1225.0 

.2554 

14.8612 

.50 

4  CO 

297.2240 

1600.0 

.3051 

12.6667 

.50 

40.0 

253.3340 

1600.0 

.3571 

10.2022 

.50 

40.0 

204.0440 

1600.0 

.4116 

7.5539 

.50 

40.0 

151.0780 

1600.0 

.4688 

4.7470 

.50 

40.0 

94.940 

1600.0 

.5288 

1.8495 

.50 

40.0 

36.990 

1600.0 

.5919 

-  1.1205 

.50 

4C.0 

-  22.410 

1600.0 

.6582 

-  4.  163o 

.50 

40.0 

-  83.260 

1600.0 

.7281 

-7.2599 

.50 

40.0 

-145.1980 

1600.0 

.2946 

15.1715 

.50 

45.0 

341  .3588 

2025.0 

•  3445 

12.7526 

.50 

4  5.0 

286.9335 

2025.0 

.3967 

10.1319 

.50 

45.0 

227.9678 

2025.0 

.4513 

7.3707 

.50 

45.0 

165.8408 

2025.0 

.5086 

4.5188 

.50 

45.0 

101.6730 

2025.0 

.5686 

1.5944 

.50 

45.0 

35.8740 

2025.0 

.6316 

-  1.4093 

.50 

45.0 

-  31.7093 

2025.0 

.6978 

-  4.4674 

.50 

45.0 

-100.5165 

2025.0 

.7674 

-  7.5479 

.50 

45.0 

-169.8278 

2025.0 

.3328 

15.2394 

.50 

50.0 

380.9850 

2500.0 

.3829 

12.6574 

.50 

50.0 

316.4350 

2500.0 

.4352 

9.9418 

.50 

50.0 

248.5450 

2500.0 

.4899 

7.1355 

.50 

50.0 

178.3875 

2500.0 

.5472 

4.2249 

.50 

50.0 

105.6225 

2500.0 

.6072 

1.2669 

.50 

50.0 

31.6725 

2500.0 

.6701 

-  1.7206 

.  50 

50.0 

-  43.0150 

2500.0 

.7361 

-  4.7692 

.50 

50.0 

-119.230 

2500.0 

.8055 

-  7.8473 

.50 

50.0 

-196.1825 

2500.0 

Table  3    (cont.). 
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THEORETICAL 

LOGARITHMIC 

in 

T* 

Tj 

^f  . 

-1.5362 

2=9957 

-     .4574 

3.5553 

-1.3295 

2.9957 

-  .5452 

3.5553 

-1.1507 

2.9957 

-  .6860 

3.5553 

-  .9924 

2.9957 

-  .8884 

3.5553 

-  .8491 

2.9957 

-1.1580 

3.5553 

-  .7181 

2.9957 

-1.5141 

3.5553 

-  .5964 

2.9957 

-1.9988 

3.5553 

-1.3649 

2.9957 

-  .3886 

3.6889 

-1.1871 

2.9957 

-  .5093 

3.6889 

-1.0297 

2.9957 

-  .6775 

3.6889 

-  .8877 

2.9957 

-  .9051 

3.6889 

-  .7576 

2.9957 

-1. 1999 

3.6889 

-  .6371 

2.9957 

-1.5882 

3.6889 

-  .5244 

2.9957 

-2.1194 

3.6889 

-1.2221 

2.9957 

-  .3466 

3.8067 

-1.0657 

2.9957 

-  .4931 

3.8067 

-  .9246 

2.9957 

-  .6820 

3.8067 

-  .7956 

2.9957 

-  .9301 

3.8067 

-  .6761 

2.9957 

-1.2514 

3.8067 

-  .5646 

2.9957 

-1.6734 

3.8067 

-  .4595 

2.9957 

-2.2557 

3.8067 

-1.1002 

2.9957 

-  .3230 

3.9120 

-  .9599 

2.9957 

-  .4873 

3.912U 

-  .8319 

2.9957 

-  .6982 

3.9120 

-  .7136 

2.9957 

-  .9676 

3.9120 

-  .6029 

2.9957 

-1.3138 

3.9120 

-  .4989 

2.9957 

-1.7702 

3.9120 

-  .4003 

2.9957 

-2.4102 

3.9120 

Table   3    (concl.). 


87 


THEORETICAL 

L< 

2GARITHMIC 

-ir, 

ir2 

Tj 

"35 

-1.5362 

3=1824 

- 

.6931 

3.5553 

-1.3295 

3.1060 

- 

.6931 

3.5553 

-1.1507 

3.0012 

- 

.6931 

3.5553 

-  .9924 

2.8679 

- 

.6931 

3.5553 

-  .8491 

2.70 

- 

.6931 

3.5553 

-  .7181 

2.4893 

- 

.6931 

3.5553 

-  .5964 

2.2093 

- 

.6931 

3.5553 

-  .4826 

1.8032 

- 

.6931 

3.5553 

-  .3750 

1.0922 

- 

.6931 

3.5553 

-1.3649 

3.2133 

- 

.6931 

3.6889 

-1.1871 

3. 1209 

- 

.6931 

3.6889 

-1.0297 

3.0058 

- 

.6931 

3.6889 

-  .8877 

2.8653 

- 

.6931 

3.6889 

-  .7576 

2.6910 

- 

.6931 

3.6889 

-  .6371 

2.4723 

- 

.6931 

3.6889 

-  .5244 

2.  1837 

- 

.6931 

3.6889 

-  .4182 

1.7642 

- 

.6931 

3.6889 

-  .3173 

1.0080 

- 

.6931 

3.6889 

-1.2221 

3.2257 

- 

.693] 

3.8067 

-1.C657 

3.1247 

- 

.6931 

3.8067 

-  .9246 

3.0023 

- 

.6931 

3.8067 

-  .7956 

2.8548 

- 

.6931 

3.8  06  7 

-  .6761 

2.6754 

- 

.6931 

3.8067 

-  .5646 

2.4505 

- 

.6931 

3.8067 

-  .4595 

2.1507 

- 

.6931 

3.8067 

-  .3598 

1.7107 

- 

.6931 

3.8067 

-  .2647 

.8969 

- 

.6931 

3.8067 

-1.1002 

3.2284 

- 

.6931 

3.9120 

-  .9600 

2.5382 

- 

.6931 

3.9120 

-  .8319 

2.9928 

- 

.6931 

3.9120 

-  .7136 

2.8412 

- 

.6931 

3.9120 

-  .6029 

2.6550 

- 

.6931 

3.9120 

-  .4989 

2.4219 

- 

.6931 

3.9120 

-  .4003 

2.1138 

- 

.6931 

3.9120 

-  .3064 

1.6546 

- 

.6931 

3.9120 

-  .2163 

.7667 

— 

.6931 

3.9120 
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APPENDIX  D 


Condensed  Laboratory  Test  Data 
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Table   4.      Condensed    laboratory  data  and   program  data   for 
multiple   regression   by    IBM   1410  computer. 


labcrat: 

:ry  pclync 

MI  AL 

r, 

Tu 

T3 

^2^3 

~7jZ 

.9520 

10.0 

.5190 

5.190 

.2694 

.9430 

10.0 

.5190 

5.560 

.3091 

.9290 

10.0 

.6560 

6.560 

.4303 

.9270 

10.0 

.5510 

5.510 

.3036 

.8890 

10.0 

.7810 

7.810 

.6099 

.8760 

10.0 

.7860 

7.860 

.6178 

.8740 

10.0 

1.040 

10.40 

1.0816 

.8880 

9.40 

.4990 

4.690 

.2490 

.8740 

10.0 

.4990 

4.990 

.2490 

.8560 

11.30 

.5010 

5.6613 

.2510 

.8470 

12.10 

.4990 

6.0379 

.2490 

.7810 

14.30 

.50 

7.150 

.250 

.7810 

15.90 

.5090 

8.0931 

.2591 

.780 

15.0 

.50 

7.50 

.250 

.7340 

17.40 

.5020 

8.7348 

.2520 

LABORATORY    LOGARITHMIC 


.0492 

2.3026 

-  .6559 

.0587 

2.3026 

-  .5870 

.0736 

2.3026 

-  .4216 

.0758 

2.3026 

-  .5960 

.1177 

2.3026 

-  .2472 

.1324 

2.3026 

-  .2408 

.1347 

2.3026 

-  .0392 

.1188 

2.2407 

-  .6951 

.1347 

2.3026 

-  .6951 

.1555 

2.4248 

-  .6911 

.1661 

2.4932 

-  .6951 

.2472 

2.6603 

-  .6931 

.2472 

2.7663 

-  .6753 

.2485 

2.7081 

-  .6931 

.3C92 

2.8565 

-  .6892 
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APPENDIX  E 

IBM  1620  Computer  FORGO  Programs 
Used  For  Theoretical  Dimensional  Analysis 
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FORGO  PROGRAM  FOR 
T.O  SOLVE  FOR  PI3  H! 


IBM  1620  COMPUTER 
:LD1NG  PI2  CONSTANT 


19 

20 

102 

21 

71 

70 

61 

60 
63 

62 
65 

64 
68 

67 


TWI 
DEL 
P4  = 
P2  = 
A  =  . 
T  W  I 
DC  = 

wc= 
oc= 

OCT 

TWO 

TPC 

Rl  = 

R2  = 

R3  = 

DTL 

Hl  = 

Yl  = 

IF( 

°RI 

,-OR 

CON 

IF( 

TPC 

GO 

CON 

IF( 

TPC 

GO 

IF( 

TPC 

GO 

IF( 

TPC 

GO 

IF( 

TPC 

GO 

CON 

DEL 

FTC 

CTC 

DHC 

CTH 

VM  = 


MO  WRY 
C  =  80. 
TC  =  2. 

TWIC/DELTC 
10, 
0273 

H=(P2) *(DELTC)+TWIC 
15. 
8. 

WC*DELTC 

=(QC)/((12.)*<3.413)) 
C=TWIC-DELTC 
=TWOC-.l 
(WC)**.4997 
(TWIC-TPC)/(TW0C-TPC) 
L0G(R2  ) 

NC=(DELTC)/(R3  ) 
(QC)/( A*DTLNC) 

(  (  .0383  )*(  TPC+DTLNO+9.  1325  )*(R1  ) 
SENSE  SWITCH  1  )20,21 
NT  102,TPC»H1»Y1 
MAT(3F10.4) 
TINUE 

(ABS(Hl-Yl) ) -5 C ) 7  0,70,71 
=TPC-1. 
TO  9 
TINUE 

(ABS(H1-Y1 ))-2  0.) 60,60,61 
=TPC-.5 
TO  9 

(ABMH1-Y1)  )-3.)62,62,63 
=TPC-.l 
TO  9 

(ABS(H1-Y1 ) )-l.)64,64,65 
=TPC-.-5 
TO  9 

(ABS(Hl-Yl) )-. 5)67,67, 68 
=TPC-.01 
TO  9 
TINUE 

CG=(QC  )/(  17.46) 
=TPC-DELCG 

={ (5. ) /(9. ) )*(FTC-32. ) 
=  (1.13+( .00678*CTC)-QCT) / ( .0176  ) 
=DHC+CTC 
•526*DHC+43.4 
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PW=( (VM+.017*DC*CTC)*( 12.*DC) )/( 10  00. ) 

PWB=(PW)*{ 3.413) 

FTH=( (9.) /(5. ) )*(CTH)+32. 

QH=QC+PWB 
DELHG=(QH)/( 17.46) 
TPH=FTH-DELHG 
WH=.l 

5  R4= (WH)**.4997 
DELTH=(QH)/(WH) 
TWCH=TWIH+DELTH 
IF(TW0H-TPH)6»4,4 

4  WH=WH+1. 
GO  TO  5 

6  CONTINUE 

R5  = (TWIH-TPH)/(TWCH-TPH) 
R6=L0G(R5 ) 

H2=(  (WH)*(R6 ) )/( .0273) 

Y2=(  ( .0383 )*( TPH- ( DEL TH)/(R6)  )+9.132  5)#R4 
18  IF(SENSE  SWITCH  2)15*10 

15  PRINT  101,WH,H2,Y2 
101  FORMAT(3F10.4) 

10  CONTINUE 

IF( (ABS(H2-Y2) )-5  0. ) 52,52,53 
53  WH=WH+2. 

GO  TO  5 
52  IF( (ABS(H2-Y2) )-20. ) 22,22,23 
2  3  WH=WH+1. 

GO  TO  5 
22  IF( (ABS(H2-Y2) )-3. )27,27,29 
29  WH=WH+.5 

GO  TO  5 

27  IF( (ABSCH2-Y2) )-l.) 16,16,17 
17  WH=WH+.l 

GO  TO  5 

16  IF( (ABS(H2-Y2) ) -. 5 ) 24 ,24 ,25 
25  WH=WH+.05 

GO  TO  5 
24  CONTINUE 

P]=(OC)/(PWB ) 

P2=(TWIH-TWIC)/(DELTC) 

P3= (WC) /( WH) 

PUNCH  28,DHC,FTC,PW,DELTC,DELTH 

PUNCH  28, PI ,P2,P3,TWIC,TWIH 

PUNCH  28,WH,WC,QC,QH,PWB 

PUNCH  28,TPC,TPH,FTH,P4,DC 

28  F0RMAT(5F10.4) 
WC=WC+1. 

GO  TO  3 
END 


93 


FCRGC  PROGRAM  FOR  IBM  1620  COMPUTER 
TO  SOLVE  FOR  PI2  HOLDING  PI  3  CONSTANT 


19 

20 

102 

21 

71 

70 

61 

60 
63 

62 
65 

64 
67 

66 


MOWRY 
C  =  80. 
TC  =  2. 

TWIC/DELTC 
0273 
15. 
.5 


TWI 

DEL 

P4  = 

A  =  . 

DC  = 

P3  = 

WC  = 

QC  = 

OCT 

TWO 

TPC 

Rl  = 

R2  = 

R3  = 

DTL 

Hl  = 

Yl  = 

IF( 

PRI 

FOR 

CON 

IF( 

TPC 

GO  TO  9 

CON 


WC*D 
=  (0C 
C  =  TW 
=  TWO 
(WC) 
(  TWI 
LOG( 
NC=( 
(QC) 
(  (  .0 
SENS 
NT  1 
MAT( 
TINU 
(ABS 
=  TPC 


IF( 

TPC 
GO 
IF< 
TPC 


TINU 
(ABS 
=  TPC 
TO  9 
(ABS 
=  TPC 


GO  TO  9 


ELTC 

)  /(  (  12.  )*(3.413)  ) 

IC-DELTC 

C-.l 

**.4997 

C-TPC)/ (TWOC-TPC) 

R2) 

DELTC)/(R3 ) 

/( A*DTLNC) 

3  83)*(TPC+DTLNC)+9.1325)*(R1) 

E  SWITCH  1)20.21 

02»TPC»H1»Y1 

3F10.4) 

E 

(Hl-Yl) )-50. )70,70,71 

-1. 


(Hl-Yl ) )-20. ) 60,60,61 
-.5 

(Hl-Yl) )-3. )62,62,63 
-.1 

(Hl-Yl) )-l. )64,64,65 
-.05 


(ABS 

=  TPC 

TO  9 

(ABS (Hl-Yl ))-.5)66,66,67 

=TPC-.Ol 


IF( 

TPC 

GO 

IF( 

TPC 

GO  TO  9 

CONTINUE 

DELCG=(QC)/(  17.4.6) 

FTC=TPC-DELCG 

CTC=( (5. ) /(9. ) )*(FTC-32. ) 

DHC=( 1.13+( ,00678*CTC)-QCT)/(  .0176) 

CTH=DHC+CTC 

VM=.526*DHC+43.4 
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15 

101 

10 

53 

52 
23 

22 
29 

27 
25 

24 


PW= ( (VM 

PW3=(PW 
FTH=( (9 
QH=QC+P 
DELHG»< 
TPH=FTH 
WH=  (WO 
R4=(WH) 
DELTH=( 
TWIH=TP 
TWCH=TW 
IF(TWCH 
TWIH=TW 
GO  TO  5 
CONTINU 
R5=(TWI 
R6=LCG( 
H2=( (WH 
Y2=( ( .0 
IF(SENS 
PRINT  1 
FCRMAT( 
CONTINU 
IF( (ABS 
TWIH=TW 


+  . 

)* 

•  ) 

WB 

QM 

-D 

/( 

•** 

QH 

H 

IH 

-T 

IH 


017*DC*CTC)*( 12.*DO )/ (10  00. ) 

(3.413) 

/(5.  )  )*(OH)+32. 

J /( 17.46) 

ELHG 
P3) 
.4997 
)/(WH) 

+DELTH 

PH)6,4,4 

-1. 


E 

H-T 

R5  ) 

)*.{ 

383 

E  S 

01, 

3F1 

E 

(H2 

IH- 


PH) /(TWCH-TPH) 

R6  )  )/ (  .0273) 

)*(TPH-(DELTH)/(R6) )+9.1325)*R4 

WITCH  2) 15*10 

TWIH,H2»Y2 

0.4) 

-Y2) )-20. ) 52,52*53 

1. 


28 


GC  TO  5 
IF( (ABS 
TWIH=TW 
GC  TC  5 
IF( (ABS 
TWIH=TW 
GC  TC  5 
IF( (ABS 
TWIH=TW 
GC  TC  5 
CCNTINU 
P1=(QC) 
P2=(TWI 
PUNCH  2 
PUNCH  2 
PUNCH  2 
PUNCH  2 
FCRMAT( 
WC=WC+1 
GC  TC  3 
END 


(H2-Y2) )-3. )22,22.23 

IH-.5 

(H2-Y2) )-l.)27»27»29 

IH-.l 

(H2-Y2) )-. 5)24, 24*25 

IH-.05 

E 

/(PWB) 

H-TWIO/  (DELTC) 

8,DHC,FTC»PW,DELTC»DELTH 

8,P1,P2.P3*TWIC,TWIH 

8,WH»WC,QC»QH»PWB 

8,TPC,TPH»FTH»P4,DC 

5F10.4) 
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APPENDIX   F 


Important  Thermoelectric   Equations 
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One  of  the  most  outstanding  thermoelectric  cooling  effects  is  the 
Peltier  effect.   The  Peltier  effect  is  noted  when  two  dissimilar  electrical 
conducting  materials  are  joined  and  a  direct  current  is  passed  through  them. 
The  Peltier  effect  is  noted  by  the  fact  that  a  temperature  difference  is 
generated  between  the  two  junctions.   The  defining  function  that  describes 
this  operation  is: 

q  =(of   -Of.  )  T  I  -  %I2R  -  K(T,  -  T  ) 
c      p    h   c  he 

where: 

Q  =  net  rate  of  heat  pumped,  in  watts  (from  the  cold  junction) 

OC  =  Seebeck  coefficients  of  the  two  materials  in  volts/  C. 

I  =  direct  current  in  amperes 

R  =  total  resistance  of  the  two  materials,  in  ohms. 

K  =  thermal  conductance  of  the  two  materials,  in  watts/  C. 

T,  =  hot  junction  temperature,  in  Kelvin 

T  =  cold  junction  temperature,  in   Kelvin 

The  first  term  is  the  Peltier  cooling  quantity,  the  second  is  the  Joule 
heat  generated  (one  half  goes  to  the  hot  junction  and  one  half  to  the  cold 
junction)  and  the  last  term  is  the  heat  transferred  by  conduction  through 
the  material. 

The  heat  liberated  at  the  hot  junction  may  be  expressed  by  the  follow- 
ing equation: 

Q  =  (c*   -CX  )T,  1  +  %I2R  -  K(T,  -  T  ) 
h      p    h  h  he 

Here  the  first  term  is  the  Peltier  heat  released  at  the  hot  junction, 
the  second  term  is  one  half  of  the  Joule  heat  and  the  third  is  the  heat 
conducted  away  from  the  hot  junction  through  the  material. 
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A  relationship  has  been  evolved  between  the  Seebeck  coefficient,  ther- 
mal conductivity,  and  the  electrical  resistivity.   This  relationship  is 
expressed  as  a  factor  Z  or  figure  of  merit. 

where: 

S  =  Seebeck  coefficient  in  volts/  C. 
P  =  Electrical  resistivity  in  ohm  •  cm. 
k  =  Thermal  conductivity  in  watts/cm.*  C. 

From  this  expression  the  dimension  of  Z  is   C 
Important  Formulas: 


AT      =  i;ZT  2 
max.       c 


Q.       =  K(AT     -  AT) 
^1  max .       max . 


AT    -AT 
max. 

max.      Z  T,  T 
h  c 


(c*   )2 

AB 


where: 

AT     =  maximum  temperature  difference  that  can  be  produced  by  a 

couple,   C. 
AT     =  operating  temperature  difference  (T,  -  T  )   C. 
k     =  thermal  conductivity  -  watts/cm*  C. 
T,  and  T  =  Hot  and  cold  junction  temperatures  -  Kelvin 
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P  =  Electrical  resistivity  -  ohm-cm. 

Q.  =  maximum  cooling  effect  -  watts. 
1  max .  ° 

0  =  maximum  coefficient  of  performance  for  maximum  cooling  - 
max . 

dimensionless. 

K  =  thermal  conductance  -  watts/  C. 

Z  =  Figure  of  merit 

°£  =   Seebeck  coefficient  of  the  materials  in  volts/  C. 
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APPENDIX  G 


List  of  Symbols  and  Explanation 


LIST  OF   SYMBOLS   USED   IN  THEORETICAL  ANALYSIS 
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Symbol 


Description 


Units 


A 
C 

h 
I 
k 
L 
m 


Nu 


Pr 


Re 


w 


"wb 


cf 


cc 


'hf 


'he 


pc 


•ph 


area 

thermal  conductance 

heat  transfer  coefficient  for  flat  plate 

electric  current  (direct  current) 

thermal  conductivity 

length 

mass  flow  rate 

Nusselt  number 

Prandtl  number 

Reynolds  number 

electric  power 

electric  power 

rate  of  heat  transfer 

cooling  capacity  per  couple 

cooling  capacity  of  module 

heat  rejected  by  module 

thermal  resistance 

cold  junction  temperature 

cold  junction  temperature 

hot  junction  temperature 

hot  junction  temperature 

plate  temperature  on  cold  side 

plate  temperature  on  hot  side 


ft. 

Btu./ft.2°F.hr. 


Btu./hr.ft.  F. 

amperes 

Btu./hr.ft.°F. 

ft. 

lbs./hr. 

dimensionless 

dimensionless 

dimensionless 

watts 

Btu./hr. 

Btu./hr. 

watts 

Btu./hr. 

Btu./hr. 

ft.2hr.°F./Btu, 

°F. 

°c- 

°F. 
°C 
°F. 
°F. 
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wa 

Twih 

Twoh 

Twic 

Twoc 

v 

V 

V 
m 

wh 

W 


At 


In 


At 


At, 


At 


he 


At 
eg 

V 

e 

Tl 
7Tk 


average  water  temperature 

inlet  water  temperature  on  hot  side 

outlet  water  temperature  on  hot  side 

inlet  water  temperature  on  cold  side 

outlet  water  temperature  on  cold  aide 

velocity 

volume  rate  of  flow 

potential  drop  per  couple 

water  flow  through  hot  junction  sink 

water  flow  through  cold  junction  sink 

water  flow  rate  in  heat  transfer  equation 

heat  transfer  coefficient  for  flat  plate 
related  to  water  flow  and  water 
temperature 

log  mean  temperature  difference  between 
plate  and  liquid 

temperature  drop  of  liquid  due  to  cooling 

temperature  rise  of  liquid 

temperature  difference  between  hot 
and  cold  junctions 

temperature  difference  across  copper 
plate  and  nylon/grease  separator 

kinematic  viscosity 

density 


Q  /P  , 
c   wb 

(Twih-Twic)/4t 


F. 


'P. 


'F. 


F. 


'P. 


ft. 

/hr. 

ft. 

3. 
/sec. 

millivolts 

lbs 

./hr. 

lbs 

./hr. 

lbs 

./hr. 

Btu 

./hr.ft, 

2oF 

w  /w, 

c  h 
Twic /At 


F. 

°F. 
°F. 
°C 

°F. 

ft.2/hr. 

lb.  /ft.3 
ra 

dimensionless 
dimensionless 
dimensionless 
dimensionless 
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Table 

5.      Bas 

ic    laboratory   test  data. 

Millivo 

Its 

Minutes    for 

300  pram  sample 

W              Wu 
c               h 

Direct   current 

Test   No. 

Twic 

At 
c 

Twih 

Ath 

voltage 

T2 

constant 

1. 

1.276 

0.046 

1.742 

0.071 

2.10 

1.64 

0.830 

1.275 

0.046 

1.742 

0.072 

2.10 

1.65 

0.830 

1.275 

0.046 

1.741 

0.073 

2.11 

1.64 

0.830 

1.276 

0.047 

1.742 

0.072 

2.10 

1.64 

0.830 

2. 

1.280 

0.046 

1.742 

0.048 

1.99 

1.11 

0.825 

1.274 

0.045 

1.741 

0.047 

1.99 

1.11 

0.825 

1.277 

0.046 

1.741 

0.048 

1.98 

1.12 

0.825 

1.276 

0.046 

1.742 

0.048 

1.99 

1.11 

0.825 

3. 

1.046 

0.047 

1.505 

0.096 

2.20 

2.29 

0.805 

1.046 

0.046 

1.508 

0.095 

2.20 

2.29 

0.805 

1.047 

0.046 

1.506 

0.096 

2.21 

2.28 

0.805 

1.046 

0.047 

1.506 

0.096 

2.20 

2.29 

0.805 

4. 

1.048 

0.046 

1.505 

0.058 

2.10 

1.38 

0.795 

1.047 

0.047 

1.508 

0.059 

2.11 

1.39 

0.795 

1.046 

0.047 

1.506 

0.059 

2.10 

1.38 

0.795 

1.047 

0.047 

1.508 

0.060 

2.10 

1.38 

0.795 

5. 

0.822 

0.046 

1.276 

0.073 

2.26 

1.78 

0.785 

0.824 

0.046 

1.277 

0.072 

2.25 

1.77 

0.785 

0.823 

0.046 

1.276 

0.073 

2.25 

1.77 

0.785 

0.823 

0.045 

1.277 

0.072 

2.26 

1.77 

0.785 

6. 

0.822 

0.046 

1.276 

0.049 

2.15 

1.18 

0.780 

0.820 

0.046 

1.276 

0.051 

2.16 

1.18 

0.780 

0.823 

0.046 

1.277 

0.051 

2.15 

1.18 

0.780 

0.823 

0.045 

1.276 

0.050 

2.15 

1.18 

0.780 

7. 

0.823 

0.044 

1.276 

0.046 

2.09 

1.09 

0.775 

0.821 

0.046 

1.275 

0.046 

2.10 

1.09 

0.775 

0.821 

0.044 

1.277 

0.047 

2.10 

1.09 

0.775 

0.822 

0.046 

1.276 

0.046 

2.10 

1.09 

0.775 

Table  5  (cont.). 
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Test  No. 


Twic 


Millivolts 


At 


Twih 


At, 


Minutes  for 
300  gram  sample 


W 


W, 


Direct  current 
voltage 


8. 


9. 


10. 


11. 


12, 


13. 


14. 


15. 


77~o  constant 


0.823 

0.046 

1.508 

0.050 

2.48 

1.24 

0.800 

0.824 

0.047 

1.507 

0.049 

2.48 

1.24 

0.800 

0.823 

0.046 

1.507 

0.049 

2.48 

1.24 

0.800 

0.823 

0.046 

1.508 

0.050 

2.48 

1.24 

0.800 

0.824 

0.046 

1.336 

0.048 

2.32 

1.17 

0.775 

0.823 

0.046 

1.336 

0.047 

2.34 

1.17 

0.775 

0.823 

0.046 

1.335 

0.048 

2.33 

1.17 

0.775 

0.823 

0.046 

1.336 

0.048 

2.33 

1.17 

0.775 

0.823 

0.047 

1.249 

0.047 

2.27 

1.13 

0.770 

0.824 

0.046 

1.248 

0.046 

2.26 

1.14 

0.770 

0.822 

0.046 

1.249 

0.046 

2.27 

1.13 

0.770 

0.824 

0.047 

1.248 

0.046 

2.27 

1.13 

0.770 

1.048 

0.047 

1.856 

0.053 

2.50 

1.25 

0.845 

1.048 

0.047 

1.854 

0.052 

2.50 

1.25 

0.845 

1.047 

0.048 

1.855 

0.052 

2.51 

1.25 

0.845 

1.048 

0.047 

1.856 

0.053 

2.51 

1.25 

0.845 

1.047 

0.046 

1.508 

0.046 

2.20 

1.10 

0.805 

1.046 

0.047 

1.509 

0.047 

2.20 

1.10 

0.805 

1.048 

0.046 

1.507 

0.046 

2.21 

1.10 

0.805 

1.047 

0.046 

1.508 

0.046 

2.20 

1.10 

0.805 

1.046 

0.047 

1.710 

0.048 

2.42 

1.21 

0.820 

1.045 

0.046 

1.711 

0.048 

2.41 

1.21 

0.820 

1.045 

0.046 

1.711 

0.049 

2.42 

1.21 

0.820 

1.046 

0.047 

1.710 

0.048 

2.42 

1.21 

0.820 

1.277 

0.047 

2.018 

0.052 

2.32 

1.18 

0.855 

1.276 

0.046 

2.019 

0.051 

2.31 

1.18 

0.855 

1.277 

0.046 

2.018 

0.052 

2.32 

1.18 

0.855 

1.277 

0.047 

2.018 

0.052 

2.32 

1.18 

0.855 

1.277 

0.046 

1.838 

0.048 

2.20 

1.10 

0.830 

1.276 

0.046 

1.837 

0.050 

2.21 

1.10 

0.830 

1.277 

0.047 

1.837 

0.050 

2.20 

1.10 

0.830 

1.276 

0.047 

1.838 

0.049 

2.20 

1.10 

0.830 
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The  past  decade  has  seen  the  development  of  new  serai-conductor  materials 
which  exhibit  a  greater  thermoelectric  effect  than  previously  was  available. 
As  a  result  the  application  of  thermoelectric  cooling  devices  may  become 
more  feasible  within  the  next  decade  or  so.   Devices  which  can  cool  and  heat 
without  moving  parts  are  already  being  manufactured. 

The  application  of  thermoelectric  cooling  may  be  feasible  in  agriculture 
for  such  uses  as  milk  pasturization  and/or  milk  cooling  coupled  with  water  or 
air  heating  or  cooling. 

Design  of  thermoelectric  devices  usually  proceeds  by  selecting  some 
cold  junction  temperature  to  be  held  constant.   Such  may  not  be  the  case 
with  some  applications  as  mentioned  in  this  study.   The  development  of  sim- 
plified engineering  design  equations  can  be  of  help  in  predicting  the  operat- 
ing characteristics  of  such  devices,  especially  if  the  operating  conditions 
are  not  constant.   It  was  decided  to  investigate  the  application  of  the 
principles  of  dimensional  analysis  to  determine  the  relationship  of  the 
important  variables  involved  in  a  thermoelectric  heating /coo ling  device. 
It  was  first  necessary  to  explore  theoretically  the  dimensionless  Pi  groups 
that  could  be  formed  in  order  to  determine  the  proper  combination  that  would 
yield  an  explicit  relationship.   The  theoretical  analysis  was  performed  with 
an  IBM  1620  computer. 

A  simple  device  was  fabricated  to  handle  liquid  in  both  heat  exchangers. 
The  heat  exchanger  design  approached  that  of  a  flat  plate  but  with  the  liquid 
entering  the  chambers  at  right  angles.   An  empirical  equation  for  flat  plate 
heat  transfer  was  used  in  the  theoretical  analysis.  Operating  characteristics 
of  the  thermoelectric  module  were  derived  from  the  manufacturers  data. 

The  current  was  held  constant  at  15  amperes.  The  theoretical  analysis 


assumed  water  inlet  temperatures  on  the  cold  side  heat  sink  of  70,  80,  90  and 
100  F.   Flow  rates  were  assumed  between  about  10  lbs/hr  to  about  50  lbs/hr. 
The  temperature  drop  of  the  water  being  cooled  was  assumed  to  be  2  F  and 
held  constant  at  that  value. 

The  data  obtained  from  the  theoretical  analysis  with  the  IBM  1620  com- 
puter was  programed  for  multiple-regression  analysis  with  the  IBM  1410  com- 
puter.  The  multiple-regression  analysis  yielded  an  equation  which  showed  the 
mathematical  relationship  of  the  Pi  groups. 

Equipment  was  assembled  in  the  laboratory  and  analysis  similar  to  that 
used  with  the  computer  program  was  followed.   The  data  was  then  programed  for 
multiple-regression  analysis  by  the  IBM  1410  computer. 

The  results  of  the  laboratory  tests  were  in  general  agreement  with  the 
relationships  as  predicted  by  the  theoretical  analysis. 


